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Abstract 

In the early stage of ultrarelativistic heavy-ion collisions chiral symmetry is restored temporarily. During 
this so-called chiral phase transition, the quark masses change from their constituent to their bare values. 
, This mass shift leads to the spontaneous non-perturbative creation of quark-antiquark pairs, which effectively 

£N| ■ contributes to the formation of the quark-gluon plasma. We investigate the photon production induced by 

this creation process. We provide an approach that eliminates possible unphysical contributions from the 
D ■ vacuum polarization and renders the resulting photon spectra integrable in the ultraviolet domain. The off- 

equilibrium photon numbers are of quadratic order in the perturbative coupling constants while a thermal 
production is only of quartic order. Quantitatively, we find, however, that for the most physical mass-shift 
scenarios and for photon momenta larger than 1 GeV the off-equilibrium processes contribute less photons 
than the thermal processes. 
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I. INTRODUCTION 



Ultrarelativistic heavy-ion collision experiments allow for studying strongly interacting matter 
under extreme conditions. Such experiments are currently performed at the Relativistic Heavy Ion 
Collider (RHIC) at the Brookhaven National Laboratory (BNL) and at the Large Hadron Collider 
(LHC) at the European Organization for Nuclear Research (CERN). Furthermore, they will be 
carried out at the future Facility for Antiproton and Ion Research (FAIR) at the Helmholtz Center 
for Heavy Ion Research (GSI) and at the future Nuclotron-based Ion Collider Facility (NICA) at 
the Joint Institute for Nuclear Research (JINR). One main objective of these experiments is the 
creation and exploration of the so-called quark-gluon plasma (QGP), a state of matter of deconfined 
quarks and gluons. The two basic features of the strong interaction, namely confinement [l|_and 
asymptotic freedom 0,0] predict that this state is created at high densities and temperatures 0-0]) 
which occur during ultrarelativistic heavy-ion collisions. 

The lifetime of the QGP created during a heavy-ion collision is expected to be of the order of 
up to 5 — 10 fm/c 0]. After that it transforms into a gas of hadrons. Thus, experiments cannot 
access the QGP directly, which makes the determination of the properties of this state difficult. 
Therefore, it is important to find theoretical signatures that provide a distinction between a hadron 
gas and a QGP. Furthermore, one has to identify experimental observables from which one can draw 
conclusions on theses theoretical signatures [6|-la| . One important category of these observables are 
direct photons as electromagnetic probes. As they only interact electromagnetically, their mean 
free path is much larger than the spatial extension of the QGP. Therefore, they leave the QGP 
almost undisturbed once they have been produced and thus provide direct insight into the early 
stage of the collision. 

One important aspect in this context is that the quark-gluon plasma, as it occurs in a heavy-ion 
collision, is not a static medium. It first thermalizes over a finite timescale, then expands and 
cools down before it hadronizes finally. This non-equilibrium dynamics has always been a major 
motivation for investigations in non-equilibrium quantum field theory Besides the role of 



possible memory effects during the time evolution 2114291] , it is of particular interest how the finite 
lifetime of the quark-gluon plasma itself affects the resulting photon spectra. 

The first investigations on this topic were done by Boyanovski et al. 3(], 31]. The authors first 



specified the density matrix at some initial time, to, for a thermalized quark-gluon plasma not 
containing any photon, i.e., 

Tr e~l 3H Qcn 

with (3 = 1/T. Afterwards, the authors propagated the system from the initial time, to, to a later 
time, t, under the influence of the electromagnetic interaction and determined the photon number 
at this point of time. One main result was the prediction of contributions of first-order QED 
processes, which are forbidden kinematically in thermal equilibrium. Furthermore, the photon 
spectra resulting from these processes flattened into a power law decay for photon energies oj^ > 
1.5 GeV (u)r = \k\ with k denoting the three- momentum of the photon) and thus dominated over 
higher equilibrium contributions such as gluon Compton scattering and quark pair annihilation 
with associated gluon production in that domain. 



On the other hand, the investigations in 30j, [3l| were also accompanied by the problem that 



the photon spectra decayed too slowly for being integrable in the ultraviolet (UV) domain. In this 
domain, these spectra behaved as 1/^ > which means that the total number density and the total 
energy density of the emitted photons were logarithmically and linearly divergent, respectively. 
Furthermore, the authors did not include the process where a quark-antiquark pair together with a 
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photon is spontaneously created out of the vacuum. Such a process is conceivable, as the temporal 
change of the background can provide energy for particle creation. 

An inclusion of this process in a subsequent work [33] revealed the even more serious problem 
that the contribution from the vacuum polarization, which is included in this process, is divergent 
for any given photon energy, wr. The authors argued that this contribution is unphysical and can be 
eliminated (renormalized) by rescaling the photon field operators with the vacuum wavefunction 
renormalization, y/~Z, but they did not provide a detailed calculation from which this could be 
inferred. Furthermore, one still encounters the problem with the UV behavior of the remaining 
contributions. 

Later on, the topic was also picked up by Fraga et al. 33|, 34] where the ansatz used in 3ol-33] 
was considered as doubtful, as it came along with the mentioned problems. In particular, the 
concerns raised in 33|, [34j were the following: 



In 3CH32I] the time at which the photons are observed has been kept finite. Either this 
corresponds to measuring photons that are not free asymptotic states or it corresponds to 
suddenly turning off the electromagnetic interaction at this point in time. Both cases are 
questionable. 

A system of quarks and gluons, which undergoes electromagnetic interactions, necessarily 
contains photons. Hence, taking an initial state without any photons and without the Hamil- 
tonian for electromagnetism corresponds to switching on the electromagnetic interaction at 
the initial time, which is questionable as well. It was shown in 35| that the ansatz used in 



3GH32I] is indeed equivalent to such a scenario. 



• The divergent contribution from the vacuum polarization is unphysical and thus needs to be 
renormalized. Nevertheless, the renormalization procedure presented in 32|] is not coherent 
since no derivation of the photon yield with rescaled field operators has been presented in 



3J. 



The authors of [33j, [34j did, however, not provide an alternative approach for how to handle the 
mentioned problems in a consistent manner. Solely in [35| it was indicated that the question of 
finite-lifetime effects could be addressed within the 2PI (two-particle irreducible) approach even 
though the conservation ofgauge invariance remains challengin g. Later on Boyanovsky et al. 
insisted on their approach j36l ] and objected to the arguments by 33|, 34] as follows. 



Non-equilibrium quantum field theory is an initial-value problem. This means that the 
density matrix of the considered system is first specified at some initial time and then 
propagated to a later time by the time-evolution operator. For that reason, the Hamiltonian 
is not modified by introducing a time-dependent artificial coupling as it would be the case 
for a 'switching on' and a later 'switching off' of the electromagnetic interaction. 

The quark-gluon plasma, as it occurs in a heavy-ion collision, has a lifetime of only a few 
fm/c. Therefore, taking the time to infinity is unphysical as this limit requires the inclusion 
of non-perturbative phase-transition effects on the photon production. 



The renormalization technique of |36| provides a rescaling of the photon field operators such 
that the photon number operator actually counts asymptotic photon states with amplitude 



one. 



This debate has actually been one of our motivations to study the mentioned problems. In 
the first attempt, we have modeled the finite lifetime of the (thermalized) quark-gluon plasma 
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during a heavy-ion collision by introducing time-dependent occupation numbers in the photon self- 
energy [37( ■ This ansatz allowed us to renormalize the divergent contribution from the vacuum 
polarization consistently. Furthermore, if the occupation numbers are switched on and off again to 
mimic the time evolution of a quark-gluon plasma during a heavy-ion collision, it also renders the 
resulting photon spectra integrable in the ultraviolet domain if one takes into account that both 
the creation and the hadronization of the quark-gluon plasma take place over a finite interval of 
time. The photon spectra, however, remain non-integrable in the UV domain if they are considered 
at a point of time, t, where the plasma still exists. Thus, the problem with the UV behavior is not 
under control for the general case. 

It is conceivable that these shortcomings result from a violation of the Ward-Takahashi identities 



within the model descriptions 30M32I. 1371 ] on photon production from an evolving QGP. Therefore, 



in the present work, we consider a conceptually different scenario, where the production of quark 



pairs and photons results from a change in the quark mass. In contrast to [37(, such a scenario 
has the crucial advantage that it allows for a first-principle description by introducing a Yukawa- 
like source term in the QED Lagrangian. The source term couples the fermion field to a scalar, 
time-dependent background field, 4>{t). Thereby, the fermions effectively achieve a time-dependent 
mass, which is compatible with the Ward-Takahashi identities. 

During the chiral phase transition in the very early stage of a heavy-ion collision the quark 



mass drops from its constituent value, m c , to its bare value, mj,. It has been shown in [38l. 1391] that 
this change leads to the spontaneous and non-perturbative production of quark-antiquark pairs, 
which in turn contributes effectively to the creation of a quark-gluon plasma. We investigate the 
photon emission arising from this creation process. In this context, the emitted photons do not 
only serve as a signature for the finite thermalization time of the QGP itself but, in particular, 
also for the nature of the chiral phase transition. The quarks can obtain energy by the coupling to 
the time-dependent source field. Therefore, photons can be produced in first-order QED processes, 
which would be kinematically forbidden in a static thermal equilibrium. We restrict ourselves to 
these first-order QED processes but maintain the coupling to the source field up to all orders. 



Similar investigations have been performed in 40(] on electron-positron annihilation into a single 



photon in the presence of a stron g la ser field, with the pair creation induced by a time-dependent 



electromagnetic background field 411 ]. 

In this context, there is a crucial difference to the approaches in l37| : There the photon 

numbers have been considered at finite times, t. In the present work we will show, however, that 
the photon numbers have to be extracted in the limit t —> oo, i.e. for free asymptotic states, which 
are the only observable ones because they reach the detectors. For this purpose, we specify our 
initial state at to — > — oo and introduce an adiabatic switching of the electromagnetic interaction, 
i.e., 

H E M(t) -> f £ (t)H E M(t) , with f e (t) = e- £ l*l and e > . (2) 
The photon numbers are then considered in the limit t — > oo and we let e — > at the end of 



our calculation. Hence, we pursue an in/out description as suggested in [33j, |34j. In this paper 
we shall demonstrate that this procedure eliminates possible unphysical contributions from the 
vacuum polarization and, furthermore, renders the resulting photon spectra integrable in the UV 
domain if it is taken into account that the mass change takes place over a finite time interval, r. 
In particular, we also show that keeping the exact sequence of limits, i.e., taking first t — > oo and 
then e — > 0, is indeed essential and that interchanging them leads to an inadequate definition of 
the photon number. In particular, we point out that it effectively comes along with a violation 
of the Ward-Takahashi identities. This again underlines the fact that a definition of physically 
meaningful photon numbers at finite time is problematic. 
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This paper is structured as follows. In Sec. HU we present our approach to photon production 
arising from the chiral mass shift. We calculate the photon yield up to first order in a e but maintain 
the coupling to the external background field up to all orders. For this purpose, we construct our 
interaction picture in such a way that it already incorporates the underlying interaction term of the 
external Yukawa field with the quarks in the unperturbed Hamiltonian of our system. The photon 
yield from first-order QED processes is then obtained by a standard perturbative calculation. We 
shall also prove explicitly that the coupling of the fermion field to the scalar background field does 
not violate the Ward-Takahashi identities. 

Before we turn to our subsequent investigations concerning photon production, we provide in 
Sec. IIIII an insertion on pair production arising from the chiral mass shift. There we will basically 
generalize the results of [38|, [39| on the quark and antiquark occupation numbers from asymptotic to 
finite times. One crucial result is that at any finite time, t, the total energy density of the fermionic 
sector features a logarithmic divergence. This is even the case for a mass parametrization, which 
is arbitrarily often continuously differentiable. On the other hand, the truly observable fermion 
occupation number, i.e. the one at t — > +00, is not divergent. These features might indicate that 
the resulting spectrum of photons radiated by the produced quarks inherits a divergence, since for 
the photons one sums over the whole history. We will see that this is in general not the case for 
the observable photon spectrum. 

In Sec. IIVI we present our numerical results on chiral photon production. There we compare 
again different mass parameterizations, m(t). For the case of an instantaneous mass shift, the 
loop integral entering the photon self-energy features a linear divergence, caused by the scaling 
behavior of the quark-antiquark occupation numbers with respect to the fermion momentum, p. 
When regulating this divergence by a numerical cutoff in the loop integral, the resulting photon 



spectra decay oc 1/^1 in the ultraviolet domain and thus feature the same pathology as in 30. 131 ] . 
When we turn from an instantaneous mass shift to a mass shift over a finite time interval, t, which 
corresponds to a more realistic scenario, the mentioned divergence in the loop integral is cured 
and, furthermore, the resulting photon spectra become integrable in the ultraviolet domain. So 
the logarithmic divergence in the fermionic energy density at finite times does not manifest itself 
in form of a similar pathology in the the total energy density of the photonic sector. 

In Sec. |V]we close with a summary and an outlook to future investigations. Technical details 
are relegated to six appendices. 

II. GAUGE-INVARIANT MODEL FOR CHIRAL PHOTON PRODUCTION 

Our starting point is the Hamiltonian in the Schrodinger picture, H§(t), governing the time 
evolution of the system 

H s (t) = H + H g (t) + H EM , (3a) 

Ho = J d 3 x ^s(x) (-17 • V + m c ) M?) + \j^ x (4V) + hip)) , (3b) 

H g (t) = g<P(t) J d 3 x i; s (x)4> s (x) , (3c) 

Hem = [ d 3 x j MjS (x)Ag(f) . (3d) 



Here Hq contains the kinetic part for the fermions and the photons, H g {t) the coupling of the 
fermions to the external source field, (p(t), and H^m the electromagnetic interaction between the 
fermions and the photons. The current operator, j^^six), is given by 

j'm,s(£) = e^sixh^six) , (4) 
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with e and 7^ denoting the electromagnetic coupling and the Dirac matrices, respectively. The 
subscript 'S' indicates that the operators are taken in the Schrddinger picture. In this picture, the 
time evolution of the system is described by the density matrix, ps(t), which reads 

^s(t) = ^2Pn(t)\n){n\ . (5) 

n 

Here {|n)} can be any orthonormal and complete set of state vectors, i.e., it has the properties 

I re) (re| = I and (n\m) = 5 nm , 

n 

and p n (t) denotes the probability that the system is found in the state |re) at a given time, t. The 
density matrix @ obeys the equation of motion 



id t p S (t) = H s (t),p s (t) 



This equation is formally solved by 

PS (t) 
^s(Mo) 



Us(t,t )Ps(to)uUt,t ) 



T < exp 



dt'H s (t') 



to 



(6) 

(7a) 
(7b) 



Here to is the initial time and T denotes time ordering. Ug(t,to) is the so-called time-evolution 
operator which solves the equation of motion 



id t Us(t,t ) = H s (t)U s (t,t ) . 
The expectation value of an observable characterized by the operator Os is given by 

(6s) = Tr{pg(t)6 s } 

= Tr ^ps(t )ul(t,to)6 s Us(t,to)} 
= Tr{p s (t )6 H (t)} • 
In the last step, we have introduced the operator in the Heisenberg picture, 

d H (t) = ul(t,t )6 s Us(t,t ) . 



(8) 



(9) 



Since the scalar background field, <f)(t), is assumed to be classical and only time-dependent, the 
fermions effectively obtain a time dependent mass, 



m(t) =m c + gcj)(t) 



(10) 



As we shall demonstrate below (see Eq. (|6ip ) and again in greater detail in appendix iBl this is 
compatible with the Ward-Takahashi identities. It has been shown in 38|, |39(] that the change 
of the quark mass from its constituent value, m c , to its bare value, rre^, during the chiral phase 
transition in the very early stage of a heavy-ion collision leads to the spontaneous pair creation of 
quarks and antiquarks. We now investigate the photon emission arising from this creation process. 
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We assume that our system does not contain any quarks, antiquarks, or photons initially. The 
initial density matrix, ps(to), is hence given by 



PS(*0) = \0gg) (0 q g\ ® |0 7 ) (0 7 | 



(11) 



Here \0 qq -) and |0 7 ) denote the vacuum states of the fermionic and the photonic sector, respectively. 
Since the fermion mass changes in time, it is important to point out that \0 qq -) is defined with 
regard to the initial, constituent mass, m c . The initial time, to, is chosen from the domain where 
the quark mass is still at this value, i.e., to < w ith t Q denoting the time at which the change of 
the quark mass begins. In the case of parameterizations, m(t), for which the time derivative, m(t), 
has a non-compact support, it is sufficient to ensure that g<j)(t) <C m c for t < t Q . 

As the electromagnetic coupling is small, we pursue a calculation at the first order in a e but keep 
all orders in g. For this purpose, we construct an interaction picture in a way that it incorporates 

Hg(t), 



6j(t) = u°'\t,t )d s u$(t,t ) , 

U°(t,t ) = rjexp -\J dt'H 9 (t') I 
H 9 (t) = Ho + H 9 (t) . 



(12a) 
(12b) 
(12c) 



Here we have introduced the subscript 'J' in order to distinguish our interaction picture from the 
standard one in which only the kinetic part of the Hamiltonian H 9 (t) is included in the time- 
evolution operator. This standard interaction picture, denoted by 'I', will also come into play 
later. It can be shown that in our interaction picture the operators obey the analogous equations 
of motion as in the standard one 



id t Oj(t) = |Oj(t),flo,j(<) 
where Ho t j(t) is given by 

H ,j(t) = U^(t,t )H 9 (t)U°(t,t ) . 
Furthermore, we can construct an interaction-picture time-evolution operator, Uj(t,to), as 

#j(Mo) = rjexp -if\t'H 3 (t') 
Hj(t) = U 9 '\t,t )H EM U^t,t ) . 
In analogy to the standard case, it fulfills the equations of motion, 

id t Uj(t,t ) = Uj(t,t )H 3 (t) , 

as well as the identity, 

Us{t,t ) = U${t,t )Uj(t,t ) . 
Accordingly, the Heisenberg-picture operator representation ([9]) is related to (|12aj) by 

6 u (t) = u}(t,t )dj(t)u 3 (t,to) . 

Moreover, (|12ap can be related to the standard interaction-picture representation as 

Oj(t) = ul(t,t )d 1 (t)u 1 (t,t ) . 



(13) 
(14) 

(15a) 
(15b) 

(16) 

(17) 

(18) 

(19) 
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Here we have introduced 



£/i(Mo) 

Hf(t) 

Oi(t) 
£>o(Mo) 



T < exp 



dt'H°(t') 



to 



g<t>(t) J d x i)i(x)ipi(x) , 
_-iH (t-to) 



and taken into account that 



U°(t,t ) = U (t,t )U 1 (t,t ) . 



(20a) 

(20b) 

(20c) 
(20d) 

(21) 



Before we can start with our calculations on photon production, we still have to determine the 
form of the fermion- and photon-field operators within the interaction-picture representation, 'J'. 
For this purpose, we take into account that in the standard interaction picture, T, the photon-field 
operator is given by 

1 [e^jfe, X)ai(k, A, t)e^ 3 + s"'*(k, X)a\{k, A, t)e~^ 



E 



(2vr) 3 JET, 



a\(k,X,t) = 4(M)e ia ^*-*°) . 
Here uj^ is the free photon energy given by 



, (22a) 

(22b) 
(22c) 



\k\ . 



(23) 



The creation and annihilation operators, a s (k, A) and as(k,X), fulfill the commutation relation, 

a s (k,X),al(k',X')\ = (27t) 3 <5 A a'<5 (3) ik-k' 



(24) 

with all other commutators vanishing. Since Ui(t,to) contains only fermion- field operators and 
hence commutes with both aj(k,X,t) and A^(x,t), it immediately follows from (|19p that 



aj(k,X,t) = ai(k,X,t) , 
A»(x,t) = A»(x,t) . 

The fermion-field operator, V>j(x), obeys the equation of motion, 

-i7o7 • V + 7om(t) ^j(x) , 
with m(t) given by ()10p . In a more compact form, this can be rewritten as 

[vfdp - m{t)\ ^j(x) = , 



(25a) 
(25b) 



(26) 



(27) 



which is just the Dirac equation with a time-dependent, scalar mass. Since both momentum and 
spin are conserved, it is convenient to expand ijjj( x ) in terms of positive and negative energy 
eigenfunctions with momentum, p, and spin, s, 

d 3 p r* *t 
7^X3 [ b f,s^p,8,tW + d_ fiS ipp, s ,i(a 



vjj(^) = Yl 



(28) 
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The creation and annihilation operators are constructed such that they obey the anticommutation 
relations 

{*>l,s>kr} = (2n) 3 5 rs 5W(j>-q) , (29a) 
{^kr} = (2n) 3 5 rs 5W(j>-q) , (29b) 

with all other anticommutators vanishing and that both b^ s and d-p jS annihilate the initial 
fermionic vacuum state, 



bps \0 q q) = d-p, s \0qq) = , (30a) 

<<Wl b p,s = <<Wl <$-p s = • (30b) 
The positive- and negative-energy state wavefunctions 3+1(2;) fulfill the Dirac equation (|27p . i.e., 

[vfdp - m(i)] ipp s , n (x) = (31) 

with the initial conditions 

^p,s,-\{x) -> ^ a ,t( x ) = Uc(p,s)e~^ E ?~ p ' g) for t < i' , (32a) 

^a^) ^ s4 (^) = «c(p, sK^ 2 * 4 ^ for t < to . (32b) 

Choosing the Dirac representation for 7^, i.e., 

with <T* denoting the Pauli matrices, the spinors u c (p,s) and v c (p,s) read 

w c (p>) = I . c v sp I , (33a) 

(sin y& x s \ 
c ff-p ■ ( 33b ) 

cos 99^. and sin y£ are given by 



E1+m c 



2E$ 



/ £"i — m c 

s ' m ^ = y^r- (34b) 

£y is the onshell particle energy for given momentum, p. The superscript 'c' denotes that it is 
defined with respect to the initial constituent-quark mass, m c . The Weyl spinors, Xs, can be chosen 
as any orthonormal set of two-component vectors fulfilling the completeness relation 



^2 XsXs = 1 ■ 



In order to solve (|31|) . we parametrize V , p,s,t^( x ) by 



^P,sM=\l M Bpl K*. (35) 
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which leads to the following equations of motion for ap^(t) and /3p^(t) 



idtap >n (t) = p/3p n (t) + m(t)ap M {t) , (36a) 
id t (3pn(t) = Paftntt) ~ m (. t )^,n( t ) • ( 36b ) 



The initial conditions read 



ap^(t) — ) 


• cos <p£ e" 


-lEIt 

V 


for i < t 




sin (yj^. e 4 


iElt 

V 


for t < t 


%,t(*) -» 


sin tp^ e" 


LESt 
p 


for t <t 


/%,+(*) -» 


- cos tp~ e 


-H-ESi 
p 


for t < i 



(37a) 
(37b) 
(37c) 
(37d) 

It can be shown that the wavefunction parameters, ap^(t) and (3p^i(t), fulfill the normalization 
condition 



Kn(t)\ +\Pp,n(t)\ =i , (38) 



as well as the relations 



ap A (t) = /% t (t) , (39a) 

= -«;-,t(*) • ( 39b ) 

Relations ()38p and (j39j) essentially show that i/y s ^(x) indeed form an orthonormal basis for given 
momentum, p, and spin, s, as they obviously imply 

^Lt^^tCz) = ^^(aOV^fc) = 1 > ( 40a ) 

^.t^)^*^^) = V^(aOVfc*.t(aO = • ( 40b ) 

Now we turn to the description of photon emission induced by the chiral mass shift. Since the 
quark mass is time-dependent only and our system is initially given by the vacuum state (fTT|) . it 
is spatially homogeneous. For such a system, the photon number at a given time, t, reads 

d 6 n 7 (t) 1 ^-^ . 

! , . ( 41 ) 

V //it 



= -^3y 2^ ^a{j(A;,A,t)aH(fe,A,f 

The sum runs over all physical (transverse) polarizations. Before we can continue with the eval- 
uation of (|4ip . we first have to clarify under which circumstances an(k,\,t) together with its 
Hermitian conjugate actually allows for an interpretation as a single-photon operator. For this 
purpose, we recall the plane-wave decomposition (|22ap which together with ()25bj) implies 

A»(x) = U$(t,to)A$(x)U 3 (t,to) 

I ^ 1 [>(fc, A)a H (fc, A, ty ts + £ **(fc, A)4(fc, A, t)e-&*] , (42) 



E 

A 



and accordingly 

MM,*) = tf](Mo)6j(M,*)#j(Mo) ■ (43) 
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The interpretation of aj(k,X,t) from (|25a|) as a single-photon operator is evident since (|25b[) de- 
scribes a free electromagnetic field. On the other hand, (|42[) describes an interacting electromag- 
netic field. Hence, the same interpretation for (|43p is not justified in general. It is, however, 
possible in the limit t — > ±00 for free asymptotic fields. Such fields are obtained by introducing an 
adiabatic switching of the electromagnetic interaction, 

H E M(t) f £ (t)H E M(t) , with f e (t) = e- £ l*l and e > . (44) 

According to the Gell-Mann and Low theorem (42|, such a switching can also be applied to con- 
struct the eigenstates of an interacting theory out of those for a non-interacting theory. Upon the 
introduction of f e (t), the interaction-picture time- evolution operator turns into 



C/j(Mo) U!{t,t )=T\exp 







T |exp 


[-'/' 




L Jto 



dt'f e (t')H 3 (t>) 



(45) 



In order obtain physically well defined results for the photon numbers, we hence have to specify our 
initial state at to — > —00 and consider eq. (|42[) in the limit t — > 00 for free asymptotic states. Such 
an approach corresponds to an in/out description as suggested in 23, [34|. Since the introduction 
of f £ (t) per se is an artificial procedure, we have to take e — > at the end of our calculation. As a 
first step, we expand au(k,X,t) to first order in the electromagnetic coupling, e: 



au(k, A, t) = Uj(t, —oo)aj(k, A, t)Uj(t, —00) 

«oj(fc,A,t) + i / d*7 6 (t0 [Fj(f'),aj(fc,A,t) 



(46) 



ai(&, A, t) — i 



g^(fc, A) 



dt' / d d x/ e (t')j /li j(f,t')e 1 ^ 



where we have made use of relations (|25p and [ai(fc, A, i), j[i t j(x, t')] = 0. Upon insertion of ([4U 
into (|4ip . we obtain for the photon yield 



2w r 



d 6 n^(t) _ 7 ^(fc) 
' fe d 3 xd 3 fc ~ (2vr) 3 U 

with the underline denoting that 



ik(xi— x 2 ) 



(47) 



d s x 



We have introduced the photon polarization tensor 



dtf £ (t) 



-r\' 



k 2 



x=± 



for {1, 2, 3} 

otherwise 



(48) 



with rf" v = diagjl, —1,-1,-1} as well as the current-current correlator 

in< M (si,s 2 ) = (jl,j(x2)jv,j( x i)) > ( 49 ) 

which describes the photon self-energy to first order in a e . The average is taken with respect to the 
initial density matrix, p(to — > —00), which is simply given by the vacuum expression (HI]). Hence, 
we have 



H1< I1 (X 1 ,X 2 ) = (0gg\ J ] j(x 2 ) ju,j(X!) \0gg) 



(50) 
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In the interaction picture representation, 'J', the current operator @ reads 

3imA x ) = e ^j( x h^j( x ) ■ ( 51 ) 
Upon insertion of (|5ip into (|50p and performing a Wick decomposition, we obtain 

iU^(xi,x 2 ) = e 2 Tr {-f fM S^(x 1 ,x 2 huSp(x2,xi)} . (52) 

This expression corresponds to the one-loop approximation, which is depicted in Fig. [TJ The 




FIG. 1: The photon self-energy, UT^ v (xi, X2), is given by the one- loop approximation with fermion propa- 
gators dressed by the background field, 4>{t). 



propagators entering expression (|52|) are given by 

Sp{xi,x 2 ) = -i^j(xi)^j(x2)y 

/d 3 p 
^^Vy,s,t(>i) ® ^p,s,^(x 2 ) , (53a) 

Sp{x 1 ,x 2 ) = 1^(3:2)^(^1)/ 
/" d 3 

= _/ (^3"fe^( Xl ) ® $P,s,l( x 2) ■ (53b) 

Since the spatial dependence of ipp ;S ,-[i(x) is included entirely in the factor e ip ' x , it is convenient to 
formally separate it off via 

^s, n {x) = 1>'a., n ®tP* ■ (54) 

Then ^ sM (t) fulfills 

[h°d t + fpi - m(t)] = ' ( 55 ) 

With the help of ([33]) . we can rewrite ([33]) as 



5|:(xi,X2) = y — %5'>(p;t 1 ,t2)e ip -^-^ , (56a) 

S£(xi,x 2 ) = y ^5<(p,t 1 ,t 2 )e i P^-^) , (56b) 
with the propagators in mixed time-momentum representation given by 

S^(p,h,t 2 ) = -i^V^ s ,t(*i)^ s ,t(* 2 ) > ( 57a ) 

s 

S^ipMM) = iX>a^(*i)$?,a(*2) • (57b) 
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By virtue of eq. (|57p we can express the photon self-energy in mixed time-momentum representation 

by 



iH< (x u x 2 ) = I -j^JL. in<,(p,t 1 ,t 2 )e i ^^-^) , (58a) 

(58b) 



/do 
Tr {7m 5 ' < (P + £ h,t 2 )j„S > (q, t 2 ,h)} . 

Finally, this enables us to rewrite the expression for the photon yield (|47|) as 

2 ^d^ = i~J tl l-J t2 ^^^ h ^^^ , (59) 



where we have introduced 



m<(k,t l ,t 2 ) = 1 ^(k)m<Jk,t 1 ,t 2 ) 



We show in Appendix [A] that ([59]) can be written as the (space-time integrated) absolute square of 
a first-order QED transition amplitude and is thus positive (semi-) definite. Therefore, it cannot 
adopt un phy sical negative values. 

As in 37J], the photon self-energy entering (|59p is given by the one- loop approximation (|58b[) . 
The crucial difference, however, is that the underlying scenario has been addressed within a first- 
principle description. In particular, the propagators entering (|58b|) are determined by the equations 
of motion, 



{h°d tl + YPi - m(h)) Sf{p, h,t 2 ) = , (60a) 
(i 7 S t2 - f Pi + m(t 2 )) Sf(p, h,t 2 ) = . (60b) 



From (|60p it follows that our description fulfills the Ward-Takahashi identities for the photon 
self-energy, 

<9 tl in< (k,h,t 2 ) - wmf^kMM = o , (6i) 

and is hence consistent with U(l) gauge invariance. A more explicit verification of (|6ip with ()58bp 
expressed in terms of the wavefunction parameters (|35p is given in Appendix [Bj Moreover, the 
observable photon yield is extracted from (|59p for free asymptotic states by successively taking the 
limits t — > oo and then e — > after the time integrations have been performed, i.e., 

2 ^d^d4 = So (2^ J_J tl y_ oo d ^^^)^fe) in ^(fc,ti,t2)e^( tl - t2 ) . (62) 

Taking into account that (|48|) has a purely spacelike structure, it follows from ()B10p that 
illy (k, t\, t 2 ) reads in terms of wavefunction parameters 

in<(Mi,t 2 ) = 4e 2 J -j^L [c^tCtO/S^CtO^^CfejQKtCfe) 



x(px + uA / (63) 



+ ^ jt (<i)«^(*i)«^(t2)^;t(t2) 
K(*i)^ + u(*0«; + ^(i2)^t(*2) 
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with p and x denoting the absolute value of the fermion momentum, p, and the cosine of the angle 
between p and k, respectively. It follows from (|37p that (|63p reduces to the vacuum polarization if 
both time arguments, t\ and t 2 , are taken from the domain where the fermion mass is still at its 
initial value, m c , 

in^ )0 (fc,ii,i 2 ) = iII^ )0 (A;,*i - t 2 ) 



2e 2 



f^P_i px(px + u, s ) + m* \ ^ + E,)( tl -t 2 ) 

Due to the chiral mass shift, (|63p will acquire an additional non-stationary contribution, 

m^(k,t 1 ,t 2 ) =m^ (k,ti - 1 2 ) + iAU^(k,ti,t 2 ) , (65) 

depending on both time arguments separately In appendix [A] we show that the contribution from 
the vacuum polarization (|64|) vanishes when taking the successive limits t — > oo and e — > so that 
only contributions from mass-shift effects characterized by iAIL^(/c, ti,t 2 ) remain. Thereby, we also 
point out that keeping this order of limits is indeed crucial to eliminate the vacuum contribution 
and that the latter shows up again if the limits are interchanged. Moreover, we demonstrate 
in appendix O that adhering to the correct order of limits is also essential to obtain physically 
reasonable results from the mass-shift effects. Together with (|63p . expression (|62p describes photon 
production induced by the chiral mass shift at first order in a e but to all orders in g. 



III. PAIR PRODUCTION FROM DYNAMICAL MASS SHIFTS 



Before we turn to the numerical investigations on photon production arising from the chiral mass 
shift, we first provide an insertion on quark-pair production. It has been shown in 38, 3^] that the 
asymptotic quark/antiquark occupation numbers are highly sensitive to the order of differentiability 
of the considered mass parametrization, m(t). We are now going to extend these investigations to 
the time dependence of the quark and antiquark occupation numbers for different mass functions, 
m(t). We consider pair production arising from the chiral mass shift only. The starting point for 
our considerations is hence the fermionic part of the Hamiltonian in the interaction picture, J, 



d 3 x ipj(x) — ry • V + m(t) "4>j(x) . 



(66) 



To simplify the notation, the subscript 'J' is dropped from now on. With the help of (J28J) and ([35 
we can rewrite (1551 as 



m = E 



d 3 p 



{m 



St Lt, 

p,s P> & 



+ A(*)^/_ p> + A*(t)d_ p> S p> } 



(2vr) 3 

where we have introduced 

n(t) = 4>p, s ,t( x ) -if • V + m{t) Vy, s ,t( x ) 

= ~ipp, s ,l(x) -i7-V + m(t) rpp, s ,i(x) 

= p [a£ t) + /% t (t)aft t (t)] +m(t) [\a At (t)\ 2 - |/% t ,(t)| 

A(t) = 4>p :S ,t(x) -if ■ V + m(t) ipp, s ,i(x) 

= P [a^ )t (t)^(t) + f3^ t (t)a p ^(t)] + m(t) [a} A (t)a p ^(t) - ^(t)^(t)] 



(67) 



(68a) 



(68b) 
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The second equality in (|68aj) follows immediately from (|39p . Introducing 



we can rewrite (|67p in an even more compact form, 



Now the particle number density is extracted from (|70p by diagonalizing this expression with 
respect to b^ s and d_p- jS via a Bogolyubov transformation [42j ]. 

/ ^ _ Z' Z' ) _ n(±\ ( b p,s ^ 



JL..M J " WW = C(f) J ' P1) 



In order to maintain the anticommutation relations (|29p under ()7ip . the Bogolyubov coefficients 
have to satisfy the relation 

|£ P >(t)| 2 + k>(*)| 2 = l • ( 72 ) 
Furthermore, they have to fulfill the initial conditions, 

£ P >(*o) = 1 , (73a) 

ift s (*o) = . (73b) 

The Bogolyubov particle number density for given momentum, p, and spin, s, is then defined as 




= (2^sEk.( t )l s 

By virtue of (|7ip . we can obviously rewrite ()67p as 



(75) 



The adjoint matrix C^(t) denotes the inverse transformation of (|71|) and reads 
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For As{t) to be diagonal, £p iS (t) and rjp jS (t) have to be determined such that 



(77) 



are the (orthonormal) eigenvectors of A(t) for the respective eigenvalues, \±(t). These are obtained 

as 



det 



A(t) - A±(t)J 







O A ± (t) = ±v/fi 2 (t) + |A(t)| 2 
With the help of (f38|) and ([39]) . they are further evaluated to 

\±{t)=E f (t) . 
Here we have introduced the dispersion relation 



E p {t) = y/p* + m?(t) . 
In terms of the transformed operators (|71|) . the Hamiltonian (|67p then reads 

d 3 p 



(78) 



(79) 



(80) 



H(t) 



E 

s 

E 



(2tt) 3 
d 3 p 



Ep(t) [&^(%>(i) - d-p, a (t)S_ Ptt (t) 



(81) 



In the second step, H(t) has been normal ordered with respect to (|71|) in order to avoid an infinitely 
negative vacuum energy. It follows immediately from (|8ip and (|7ip that the energy density is given 
by 



d 3 x 



j^3 E p(t) | 1 2 



(82) 



We see that it corresponds to (|74p integrated over the momentum modes, p, which justifies the 
definition of (|74|) as the particle- number density at a given time, t. The additional factor of 2 
in (|82p arises from the fact that this expression describes the energy density carried by quarks 
and antiquarks together whereas (|74|) corresponds to the number of quarks which is equal to the 
number of antiquarks. Together with 



A(t) - A±(t)/J v±(t) = ±Ep(t)v ± (t) , 
we obtain the following linear system of equations for £fi s (t) and r)p )S (t), 

+A*(t)»ft s (t) = E p (t)^ tS {t) , 
A(t)fo,(*) - n(t)r] fiS (t) = Ep(t)rfr s (t) . 



(83) 



(84a) 
(84b) 



16 



In appendix [D] it is shown that this system is solved by 



e P >(*) 



e iE f (t)t 
e iE p (t)t 



IEp{t)+m(t) 
2Ep{t) 

lEp(t) + m(t) 
2%(t) 



a p,j.(*) + 



lE0(t)-m(t) 
2E p {t) 

!Ep(t)-m(t) 
2E p {t) 



/%,t(*) 



(85a) 
(85b) 



The phase factor, e lE p^*, has been introduced to satisfy the initial condition (|73a|) . Furthermore, it 
allows us to rewrite ip(x) in terms of positive- and negative-energy wavefunctions of the respective 
momentary mass, m(t), i.e., 



d 3 p 



(27T) 



V>(x 

with V>p> t 4(a0 given by 

In analogy to (|34|>. we have introduced 



(86) 



COS^(^Xa \ 

sin^(t) ^x s y 

sin^(t) 
-cosv?p-(t) 



I e +i%(*)* e ip-i 



(87a) 
(87b) 



COS </3p(*) 

sin (/?p-(t) 



'.%ft) + m(t) 
2E p (t) 

lE p (t) - m(t) 
2Ep(t) 



(88a) 
(88b) 



The Bogolyubov transformation (|74p hence corresponds to a reexpansion of the fermion-field op- 
erators in terms of the instantaneous eigenstates of the Hamilton-density operator 



ho(t) = -i7o7 • V + lom(t) 



(89) 



which is demonstrated in greater detail in appendix [Dj The same procedure has been applied in 
. The crucial difference to our approach is that the authors us an expansion in the form of (|86p 
to derive the equations of motion for the field operators (|7ip and eventually a kinetic equation for 
the Bogolyubov particle number density (|74p from the Dirac equation with a time-dependent mass. 
To the contrary, we extract the Bogolyubov parameters and hence the particle number density by 
translating (fTTj) into relations between ipp, s ,^i(%) and V'p.s,t4.( x ) an d projecting the Bogolyubov 
parameters out of the latter (also see Eqs. ()D2|) - ()D4|) in appendix [D]) . 

With the help of (|39p . both £g s (t) and r]p !S (t) can be expressed alternatively in terms of (complex 
conjugated) negative- and positive-energy wavefunction parameters, respectively. The Bogolyubov 
particle number density (f74"|) thus reads 



d% g -(t) _ 1 
d 3 xd 3 p (2vr) 



1 + 



2pRe 



«;^(*)^(*)J + m (t) [I«pu(*)| 2 - ifcM 



(90) 
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We see that it can be expressed entirely in terms of the negative-energy wavefunction parameters, 
a P,l(t) an d fip,\{t)- This is a result one would also expect intuitively since ipp iS ,i( x ) describes 
initially a negative-energy state but then also acquires a positive-energy component from the mass 
shift measured by ([90]) . This expression is, indeed, just the absolute square of the projection of 
i^p jS ,i(x) on V , p,s,t( a; ) summed over the spin index, s. For completeness we mention that (|90|) can 
also be obtained by projecting the respective propagator (|53b[) on the positive-energy wavefunction 
of the respective current mass, m(t), 



d 6 nqg(t) 

d 3 xd 3 p 



Y^p ^2 J d 3 xi J d 3 x 2 4>l S!t (x 1 )S^(x 1 ,X2ho4>p,s^(x2) 



(2tt) 3 V 



(91) 



t!=t 2 =t 



which has been used in [38|, [39J in the asymptotic limit t — > oo. Thus ([90]) generalizes the result 
therein to finite times, t. 

For our investigations on the time dependence of (|90p . we model the change of the fermion 
mass from its initial constituent value, m c , to its final bare value, m&, by three different mass 
parameterizations, 



mi(t) 
m 2 (t) 

m 3 (t) 



m c + rrif, m c — nib 



m c + nib m c 



sign(t) , 
sign(t) ( 1 — e 



-2\t\/r 



m c + mt, m c — nib ( 2t 

tanh — 



(92a) 
(92b) 
(92c) 



with sign(t) given by 



sign(t) 



1 t > 
-1 t < 
t = 



The mass parameterizations (|92ap - ()92cp are depicted in Fig. [2j Analogously to [38j, [39J, we have 
chosen m c = 0.35 GeV and nib = 0.01 GeV and assumed a transition time of r = 1.0 fm/c. 



> 
o 



0.4 ,— 

0.35 - 

0.3 - 

0.25 - 

0.2 - 

0.15 - 

0.1 - 

0.05 - 

- 



m,(t) 
m 2 (t) 
m 3 (t) 



-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 
t [fm/cl 



FIG. 2: During the chiral phase transition, the quark/ antiquark mass changes from its constituent value, 
m c , to its bare value, rrib- 
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It has been shown in [38l . |39j that the asymptotic occupation numbers are very sensitive to 
the order of differentiability of the considered mass parametrization, m(t). For the case of an 
instantaneous mass shift described by (|92a|) . the equation of motion (|3"T|) for the positive- and 
negative-energy wave function is solved analytically with the ansatz 



w*) = rH x) b for t<0, 

/ / ^ f^iW for 1 < > 

Wffs I {%) = { \ u 

P " 1 W,»4 b frJx) + 6 Rs ^(x) for t>0. 



(93a) 
(93b) 



For t < 0, ijjp iS ^(x) and ijjp tS ,i{ x ) describe positive- and negative-energy states of mass m c , respec- 
tively, whereas they turn into superpositions of positive- and negative-energy states of mass for 
t > 0. From the continuity condition 



Vy,s,nO?,o ) = ipp SjU (x,o + ) 



we obtain 



p,s 



p,s 



1p,s 

-5. 



p.s 



cos ip b p cos (fp + sin (pp s'm ipp 
sin <p h p cos (fp — cos ^sin cpp 



(94) 

(95a) 
(95b) 



As the coefficients ap >s and f3p >s do not explicitly depend on the spin, s, this index will be omitted 
from now on. The occupation numbers thus read 



dX g -(t) 
d 3 xd 3 p 



2/3| 
(2tt) 3 

1 

(2^)3 



p + mbrric 
p p 



(96) 



for t > 0, whereas they vanish for t < 0. For p 3> m&, m c , the expression (|96|) can be approximated 
as 



d 6 n q g(t) (m c 



d 3 xd 3 p ' (2tt) 3 2p 2 



+ 0{l/p A 



(97) 



which means that the total particle number density and the total energy density of the fermionic 
sector are linearly and quadratically divergent, respectively. 

This artifact can be removed if the mass shift is assumed to take place over a finite time interval, 
r. In this case, the occupation numbers are obtained by solving (|36p numerically for the negative- 
energy wavefunction parameters, ap^^t) and f3p^(t), which are then inserted into (|9Up . Fig. [3] 
compares the asym pto tic p article spectra for the different mass parameterizations, mi(t). 

Analogously to [38l . [39J, we find that if we turn from m\(t) to m,2{t), which is continuously 
differentiable once, the occupation numbers decay oc 1/p 6 and are hence suppressed relative to 
the case with the instantaneous transition. Moreover, if we turn from rri2(t) to ms(t), which is 
continuously differentiable infinitely many times, the occupation numbers are further suppressed 
to an exponential decay. In the limit r — > 0, both rri2(t) and ms(t) reproduce expression (|96|) . 
which is depicted in Fig. [H 

We shall briefly explain how the sensitivity of the asymptotic occupation numbers on the mass 
parametrization, rrii(t), comes about. For this purpose, we consider the Bogolyubov parameters in 
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FIG. 3: Asymptotic particle spectra for the different mass parametcrizations given in Eqs. (1921) . The decay 
behavior is highly sensitive to the order of differentiability of m(t). Both for rri2(t) and m^t), we have 
chosen r = 1.0 fm/c. 




FIG. 4: Asymptotic particle spectra for m.2(i) (left panel) and m^(t) (right panel) for different transition 
times, t. In each case, the suppression at large p with respect to an instantaneous mass shift is the stronger 
the more slowly (r increasing) the mass shift is assumed to take place. As expected, both parametcrizations 
reproduce expression (|M| in the limit r — > 0. 



terms of negative-energy wavefunction parameters, 



E?{t)-m(t) Ep(t)+m(t) 
W* } " V 2£^j °^ (t) " V 22^t) ^ (t) ' (98a) 



E^t)+m{t) Epjt) - m(t) Q 

= y 2JW +y 2%(t) ' (98b) 

with (|98a[) following from (|85a|) and (|39|) . We do not take into account the phase factor, e iE ^ t ' t , as 
it drops out when taking the absolute square of (|98b|) to obtain the Bogolyubov particle number 
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density (f7i|) . It follows from ([36]) that G s (£) and rj^ s (t) then obey the equations of motion, 



. pm(t) 



Vp, s (t) - E f (t)Q s (t) 



In the limit p 3> m(t), these equations of motion are approximately solved by 



(99a) 
(99b) 



' dt'm(t')e 2ipt 



(100) 



2 P 

Hence, the Bogolyubov particle number density in that domain reads 



d%,(£) (m c - m b f 



d 3 xd 3 p ' (2vr) 3 2p 2 



dt' X (t')e 



tsJiipf 



where we have formally introduced x(t) by means of the relation rh(t) = x{~t){ m c 
to — > — oo since rh(t) — > for £ — > dboo. In particular, in the asymptotic limit, £ - 
the quark/antiquark occupation numbers, 



d% 9 -(£) 



(2^) 3 2p 2 



d£'x(t')e 



(101) 

- m&) and taken 
oo, we have for 



(102) 



So for a mass shift over a finite time interval, r, (|97p is effectively modulated with the absolute 
square of the Fourier transform of %(£) from £ to 2p. So the particle numbers for p S> m c 
are suppressed by an additional factor of 1/p 2 each time the order of differentiability of m(t) is 
increased by one. In particular, for m(£) = m2(£) we have 



d Q n qq (t) 



d 6 xd A p 



t—^oo 



(2tt) 3 2 P 2 (l+p 2 T 2 ) 2 



(103) 



As a next step, we extend the investigations in [38|, [39J to the time dependence of the occupation 
numbers. Fig. [5] shows the time evolution of (|90p for different momentum modes. Here we see 
that for hard momentum modes, the occupation numbers exhibit a strong 'overshoot' over their 
asymptotic values by several orders of magnitude in the region of strong mass gradients. This means 
that the particle spectra exhibit their decay behavior characteristic for the order of differentiability 
of m(£) only in the limit £ — > oo, which is depicted in Fig. [6l 

At finite times, however, the particle spectra decay as 1/p 4 both for m,2(£) and for m^{t) for 
p ^> nib, m c . This can be understood by taking into account that at intermediate times the time 
integral in (jlOip runs from — oo to £, so that one effectively carries out a Fourier transform over 
a discontinuous function. Hence, expression (|9"U|) picks up an additional factor of 1/p 2 compared 
to the case of an instantaneous mass shift. This implies that at finite times, £, the total particle 
number density is finite whereas the total energy density features a logarithmic divergence. This is 
depicted in Fig. [7] showing its time evolution with the loop integral entering (|82p being regulated 
by a cutoff at different values of p = Ac. The considered values of Ac follow an exponentially 
increasing sequence given by 



0.2 ln^ 



A c (i) = Aie u ^ llls r 
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FIG. 5: Time dependence of occupation numbers for m2{t) (left panel) and m^i) (right panel) for r = 1.0 
fm/c. For large values of p, they exhibit an 'overshoot' over their asymptotic value around t = 0, which is 
particularly distinctive for 7713 (t). 




FIG. 6: Particle spectra for rri2(t) (left panel) and mz(t) (right panel) at different times, t, for a transition 
time of t = 1.0 fm/c. For both parametcrizations, the characteristic decay behavior only emerges for t — > oo, 
whereas one encounters a decay behavior oc l/p 4 for large p around t = 0. 



with Ai = 1 GeV, A2 = 10 GeV, and i going from to 5. This choice is such that in the 
region of strong mass gradients, the total energy density always increases by a constant amount for 
consecutive Ac(£), which reflects the logarithmic divergence of energy density (|82p . This divergence, 
however, only shows up in regions of strong mass gradients and disappears again as soon as the 
fermion mass has reached its final value. Among other things, this can also be inferred from Fig. 

El 

For m(t) = mzit) and t > 0, the particle spectrum does not decay strictly monotonously for 
p S> rrib,m c , but instead exhibits an oscillatory behavior. This can be understood by taking into 
account that 



dt' X (t': 



2ipt' 



~ 2 '/ r 1 - e~ 2t / T (cos 2pt - pr sin 2pt) 

' T\ , „9_?\9 V-W4) 



4(1+ P 2 T 2 ) (1+ P 2 T 2 ) 2 



for t > from which the oscillatory behavior of (jlOip in p at positive times becomes apparent. For 
t — > 00, the oscillating terms disappear so that the asymptotic particle spectrum shows a strictly 
monotonous decay again. 

For the sake of completeness, we also consider the case where the fermion mass first undergoes a 
change from its constituent value, m c , to its bare value, nib, and then again back to its constituent 
value, m c , after a certain period of time, tl, to simulate the temporary restoration of chiral 
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FIG. 7: Time evolution of the total energy density for m(t) = m,2(t) (left panel) and m(t) = ms(t) (right 
panel) for r = 1.0 fm/c. It features a logarithmic divergence in the regions of strong mass gradients, which 
disappears again as soon as the quark/ antiquark mass has reached its final bare value, m b . 



symmetry during a heavy-ion collision. Again we consider different mass parameterizations, 



fh\(t) 
m 2 (t) 

rh 3 (t) 



m c + m b m c -m b . / r| 2 

sign — — t 

2 2 6 V 4 

m c + m b m c -m b . / t l 

sign t H 

2 2 & V 2 



1-e 



-\2t+r L \/r 



xsign 



TL 



m c + m b m c -m b ( 2t + r L \ (r L -2t 
-tann tann 



(105a) 

(105b) 
(105c) 



Here tl denotes the lifetime of the quark-gluon plasma during which the chiral symmetry is re- 
stored. For our numerical investigations, we choose tl = 4 fm/c for which (|105ap - (|105cp are 
depicted in Fig. [SJ 



> 
6 




FIG. 8: The mass is first changed from m c to m b and then back to m c again to simulate the temporary 
restoration of chiral symmetry. 



For the case of rhi(t), where both mass shifts are assumed to take place instantaneously, the 
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Dirac equation (|3ip is solved with an ansatz similar to (J93J): 



>^, t (x) for t<-T L /2, 

a &i> h p >s ^ x ) + P^sjx) for ~T L /2<t<T L /2 

for t<-T L /2, 
7$4,8d?) + S f$},,,^ x ) for -r L /2<t<r L /2 , 



(106a) 



(106b) 



The successive application of the continuity conditions i/jp ;S ^(x, —t/2 ) = ipp, s ,^i(x, — t/2 + ) and 
'0p, s ,n(^ r / 2 ~) = ^p,s,u{x,t/2 + ) leads to 



e i(^-£»)7i/2 



-<5t 
p 

e i(i?£+M)r L /2 



2ilm{at^} 



cos </£cos y& + sin y&sin <^ 



sin ippcos (pp - cos <^>sin <pp 



Hence, for —tl/2 < t < tl/2, the occupations numbers are given by 
have 



(107a) 

(107b) 
(107c) 
(107d) 

whereas for t > tl/2 we 



d% 9 -(t) 
d 3 xd 3 p 



(2vr) 3 
2 

(2^)3 



p + mbm c 
p v 



sin 2 [E%r L 



m b ) 



(2tt)V 



sin 2 (prL) for p 3> m{„ m c 



(108) 



Consequently, for given momentum, p, the final occupation number is highly sensitive to tl and 
becomes maximal if the condition 

2<n -\- 1 

E i> t l = — 5 — 71 with n G N ° ' 

is satisfied. Comparing the asymptotic particle number density for different mass parameterizations 
rh(t), we observe the same sensitivity on the respective order of differentiability for large p as we 
did in the first scenario. This is depicted in Fig. [9] and follows immediately from expression (jlOip . 

Similarly to the first case, both (|105b|) and (|105c|) reproduce the particle spectrum of (|105a|) in 
the limit r — > 0, which is shown in figure [T0l 

Fig. [Til shows the time evolution of the particle number density for rh^{t) with r = 1 fm/c and 
tl = 8 fm/c. As it should be, the time dependence in the second scenario is the same as in the 
first one until we start changing the fermion mass back to its initial constituent value, m c . Hence, 
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FIG. 9: Asymptotic particle spectra for the different mass paramcterizations in Eqs. f|105|) . Analogously to 
the first scenario where the mass is not restored to its initial value, the decay behavior is highly sensitive to 
the order of differentiability of m(i). Both for m2{t) and for m 3 (i) we have chosen t = 1.0 fm/c. 




FIG. 10: Asymptotic particle spectra for rri2(t) (left panel) and 7713(f) (right panel) for different values 
of t. Analogously to the first scenario, the suppression at large p compared to the instantaneous case is 
the stronger the more slowly (r increasing) both mass changes are assumed to take place. As it must be, 
expression (| 108[) is reproduced in the limit r — > for both parameterizations. 



the occupation numbers first saturate at the same value as they did within the first scenario. As 
soon as a second mass gradient shows up, they are again changing by several orders of magnitude 
before they saturate at their final asymptotic value, which is usually different from the first one. 

Such a behavior follows immediately from expression (jlOip . At times when the fermion mass 
is at its bare value, nib, the integral entering it effectively represents a full Fourier transform over 
an at least continuous function while we again encounter a Fourier transform over a discontinuous 
function again when the fermion mass is being changed back to its constituent value, m c . 

For completeness, we still have to investigate in more detail how the occupation numbers are 
modified if the quark mass is changed back to its initial value, m c . For this purpose, we take into 
account that our mass parametrization (|92p and (|105|) can be written in the general form 

m(() _ rn^nv, _ ^_^ f(t) (109a) 

m = - ^i^f (« + t) ' (y - ') ■ {mb) 
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FIG. 11: Time evolution of the particle number density for rh^(t) with r = 1.0 fm/c and tl = 8.0 fm/c. They 
are again changing by several orders of magnitude when the mass is being restored to its initial constituent 
value, m c . 



with f(t) increasing monotonously from —1 to 1 and fulfilling the condition f(—t) = —f(t) (odd 
under time inversion). Hence, we have 

«'> = '(t-0*(t + ')-*(t-'Mt + 0- < 110 > 

The Fourier transform of (jllOp from t to 2p is given by 

/oo f roa \ 

dt'x(t')e 2ip *' = 2ilm e^ / dt'f (r L - 1') X (t')e 2ipt ' . (Ill) 
-oo I J — oo J 

Thus, unlike the case where the fermion mass is changed instantaneously, the occupation numbers 
for p 3> m&, m c in the second scenario are generally different from those in the first scenario 
multiplied by a factor of 4sin 2 (pr^). Nevertheless, if tl is significantly larger than r, we have 
fijL — t) ~ 1 for those times, t, where x(t) 1S significantly different from zero. We can hence 
approximate 

/oo poo 
dt'x(t')e 2ipt ' « -2isinp L r / dt' X (t')e 2ipt ' , (112) 
-oo J — oo 

where we have also taken into account that x{~ t) = x{t)- The latter follows immediately from 
f(—t) = —f{t). Hence, for r <C tl, the asymptotic occupation numbers for hard quark-antiquark 
pairs are given by 



i6 r /.^ /„_ _ \2 poo 2 



d b n q g(t) 



t—too 



(2tt)V 



dt'x(t': 



3 2ipi' 



sin 2 (pr L ) , (113) 



which just corresponds to (jlOip modified by a factor of 4 sin 2 {jpri). This is illustrated in Fig. [12] 
for the asymptotic particle spectrum of i^W, where the dotted line represents the spectrum for 
777.2 (£) multiplied by 4sin 2 (T7T£,). 

To summarize, we have found that the occupation numbers in the asymptotic limit show a 
strong dependence on the 'smoothness' of the considered mass parametrization, m(t). For the case 
of an instantaneous mass shift, they scale oc (m c — m^) 2 /p 2 for p ^> m c , mj, which means that the 
total number density is linearly divergent. This artifact is removed if the mass shift takes place over 
a finite time interval, t. In particular, the quark/ antiquark occupation numbers at p 3> ttt, c , m b and 
p> 1/t scale oc (m c — m,b) 2 /p 6 T for 777,2^) (continuously differentiable once) and are suppressed 
even further to an exponential decay for m^{t) (continuously differentiable infinitely many times). 
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FIG. 12: Comparison of asymptotic particle spectra for rri2(t) and rri2(t) for r = 1.0 fm/c. For large values 
of p, the asymptotic occupation numbers of fhi (t) essentially differ from those of m,2 (t) by a modification 
factor of 4 sin 2 ptl- 



Our investigations have shown that the pathological scaling behavior for an instantaneous 
mass shift essentially results from high-momentum Fourier components. These components are 
suppressed for mass shifts over a finite time interval such that the scaling behavior of the 
quark/antiquark occupation numbers becomes physically reasonable. A very similar effect occurs 
when considering back-to-back particle- antiparticle correlations in hi gh energy nuclear collisions 
for mass parameterizations with different order of differentiability [44J449I] . 

Furthermore, if the quark mass is changed back from its bare value, mt,, to its constituent value, 
m c , the lifetime of the chirally restored phase enters the asymptotic occupations numbers at large 
p S> mt,,m c in the form of a factor of 4sin 2 (pT^) if r <C tj_,. From the semiclassical point of view, 
the oscillations in the particle spectra emerging from this factor correspond to a multiple scattering 
of the quarks/antiquarks at both mass gradients. 

At intermediate times, however, the quark/antiquark occupation numbers decay oc 1/p 4 for large 
p in regions of strong mass gradients. This, in turn, leads to a transient logarithmic divergence in 
the total energy density of the quarks and antiquarks. At first sight such a divergence might be 
disturbing. We stress, however, that only the asymptotic energy density, i.e. for t — > +00, consti- 
tutes an observable in the sense of S-matrix theory, as the interpretation of (|74p as quark/antiquark 
occupation numbers is only justified for asymptotic times where m(t) = 0. The reason is that the 
dispersion relation (j80|) then actually characterizes free and thus detectable particles wheres it only 
describes quasiparticles for m(t) ^ fiil ]. Moreover, the asymptotic value does not show any diver- 
gence as long as the mass shift is smooth enough, which is a physically reasonable condition. Yet, 
there is one more twist in this argument when we proceed to the photon production. Photons are 
produced (and, in principle, leave the system) at any instant of time. Thus one might suspect that 
the tremendously large intermediate fermion numbers leave their imprint on the number of emitted 
photons. It is conceivable that the asymptotic number of photons, which is an accumulation over 
the whole time history, becomes large or even diverges just because there have been very many 
fermionic emitters at intermediate times. We will see in the following that this is not the case. It 
is important to understand that the disappearance of the large fermion numbers is not a damping 
effect, but a quantum mechanical interference effect. Collisions between the fermions, which would 
provide a loss rate, i.e. damping, are not included in our approach. On the other hand, the full 
quantum effects are retained, which can lead to interference patterns that are unintuitive from a 
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classical point of view. This applies in particular to quantities which are not observable anyhow 
like, for instance, the number of fermions at finite times. In the same way the asymptotic photon 
number can turn out to be reasonably small in spite of the fact that the emitted photons seem to 
pile up during the whole history of the process. 



IV. PHOTON PRODUCTION 



We now turn to our investigations on photon production. As it has been shown in the previous 
section, the asymptotic quark/ antiquark occupation numbers exhibit a strong sensitivity to the 
order of differentiability of the time dependence of the mass, m(t). In particular, they are rendered 
integrable in the ultraviolet domain, if the mass shift is assumed to take place over a finite time 
interval, r. We now investigate whether the resulting photon spectra exhibit a similar sensitivity 
and if the model of chiral photon production is accordingly suitable to describe finite-lifetime effects 
on the photon emission from a quark-gluon plasma. As one lesson from section IIIII we recall that 
only the asymptotic particle numbers constitute observables [43|, while quantities defined by the 
analogous expressions with interpolating fields have no definite interpretation as particle numbers. 
Thus, in the following we concentrate on the asymptotic photon numbers (|62p . 



A. Instantaneous mass shift 

First, we consider photon production for an instantaneous mass shift at t = as this special case 
still allows for an almost complete analytical treatment. In particular, the individual contributions 
to the photon yield allow for an interpretation as first-order QED processes and their interference 
among each other. 

On the other hand, we show that the assumption of an instantaneous mass shift comes along with 
essentially three unphysical artifacts. In section IIV A 1| we calculate the individual contributions 
to the photon yield (|62p . Thereby, we will demonstrate that the overall loop integral entering 
(|62p features a linear divergence caused by the decay behavior of the quark/antiquark occupation 
numbers ()96p for p S> m&, m c . We regulate this divergence by cutting off the loop integral at 
P = Ac- 

We will see in section HV A 21 that then the resulting asymptotic photon spectra decay as 
with cl>£ denoting the photon energy (|23|) . in the ultraviolet domain for given value of Ac*. The 
total number density and the total energy density of the emitted photons are hence logarithmically 
and linearly divergent, respectively. 

Finally, as we mention at the end of section IIV A 21 and discuss in greater detail in appendix [Fl 
another problem appears in the limit m& — >■ 0. In that limit, the loop integral over the contribu- 
tions describing quark and antiquark bremsstrahlung and quark-antiquark-pair annihilation into a 
photon feature a collinear and an anticollinear singularity, respectively. 



1. Evaluation of contributions to the photon yield 



To evaluate the contributions to the photon yield for an instantaneous mass shift, we first undo 
the contraction 

m<(k,t 1 ,t 2 ) = ^ v (k)u< fl (k,t 1 ,t 2 ) , 
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and express the photon self-energy in terms of positive and negative energy wavefunctions, i.e. 

r (ftp r _, , "IT-/ / 



r. s 



(114) 



Upon insertion of (|1 14j) into (|62p and interchanging both time integrations with the loop integral, 
we can rewrite (1621) as 
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d 3 xd 3 A; e^o (27r) 3 



d b n 
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d P re.* 



(27T) 



I £ '*(p,k,r,s)I £ u (p,k,r,s) . 
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Here we have introduced 



I*(p, k,r,s) = e I it ^(th^+^it)^ > 



(115) 



(116) 



with the underline denoting that the time integral is regulated by the convergence factor / E (t) 
e~ E l*l Moreover, with the help of (|93p and (|95ap . we obtain 



u c (p,r)7^v c (p + k,s) 



e + iwf (p, k) 



apa^u b (p, r)j^v b (p + k, s) 



1 



- apPp+tftbip, r^Ubip + k, s) 
+ Ppa^Vbip, r^Vbip + k, s) 



PpPfriFbtf, rhfjiUbip + k, s) 



e — ioj\{p, k) 
1 

1 
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1 

e + \u\(p,k) 



(117) 



To keep the notation short, we have introduced the frequencies 

b.c 

h 

u b 2 ' c (p,k) = E b f z -E b f-u 



uj b /(p,k) = E^ + E /-^ 



p+k 
b,c 

ib,c 

p+k ~P 



J/(p,k) = E b :^-E b J c + uj k , 



(118a) 
(118b) 
(118c) 
(118d) 
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Since (|118ap - (|118d|) are either positive or negative definite for both m c and mj taking the limit 
e — > leads to 

1 



lim I*(p, k, r, s) = J M (p, k, r, s) = -ie u c {p, r)j^v c {p + k, s) 



1 



u)l(p,k) 

apa^Ubip, r)j^v b (p + k, s) 

k) 

ap^ + kU b {p, r)jf,u b (p + k, s) — 
PpOt 0+ kVb(p, r)j^v b (p + k, s) 

u°{p, k) 



(119) 



l(p,k) 



With the help of 



In(p,k,r,s) = ie 



u c {p, s) = aj/u b (p, s) + /3j/v b (p, s) , 
v c {p, s) = apv b (p, s) - /%u 6 (p, s) , 

we can rewrite (|119p in the following more compact form: 

(*pa f+ tfib(p, rh„v b (p + k,s)( — - — 

+ appfrtfibfa r^Ubip + k,s)l — + — 

\u%(p,k) wi(p,fc) j 

+ PfV p+k -v b (p, r) ltl v b (p + k,s)l — - — 
+Pp a i ; + kVb{p, r)~fpU b {p + fc, s) 

Since we have ap — > 1 and 



(120a) 
(120b) 



(121) 



y uj b 4 (p,k) u{{p,k) 



Since we have — > 1 and /3 p — > for m& — > m c , expression (|121|) vanishes in this case and we 
will have no photon production, as it should be. Taking a closer look at the spinor structure 
of the particular contributions to (|121|) allows us to interpret them as first-order QED-transition 
amplitudes. It is hence convenient to split up ()12ip as 

1 



In(p,k, r,s) 


i=i 


lji(p,k, r,s) 
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(122a) 
(122b) 
(122c) 
(122d) 
(122e) 
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with the individual contributions describing the spontaneous creation of a quark-antiquark pair 
together with a photon (i = 1), quark bremsstrahlung (i = 2), antiquark bremsstrahlung (i = 3) 
and quark-antiquark pair annihilation into a photon {i = 4) . With the help of (|122p , we can rewrite 
(fTTSl) as 



d 6 n 7 

2w ■ 



7^(fc) 
' k d 3 xd 3 k ~ (2vr) 3 



I Y^I^*(p,k,r,s)P v {p,k,r,s) 



+ 2Re 



X] J M* (P' ^' r > S ) J ^ & r, s) 



i<j 



(123) 



The first term in (|123|) describes the direct contributions from first-order QED processes whereas 
the second one describes the interference among them. For further considerations, we introduce 
the shorthand notation 



hj (p, k, t) = >f (k) Re | J2 (P. r,s)4 (p, k, r, s) J , 



(124) 



in which (| 123j) reads 
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(125) 



i=i 



The evaluation of the individual contributions to (j!23|) is a lengthy but straightforward proce- 
dure and demonstrated exemplarily in appendix [Dj The direct contributions from first-order QED 
processes read 
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(126a) 
(126b) 
(126c) 
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As in section [III x denotes the cosine of the polar angle between p and k. Moreover, for the 
interference contributions we obtain 



a 2 m b\P+ k \ 



Iu(P,k)=-4e a^zP^za 



p+k H p+k V -pb 
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\wf(p,A;) u>l(p,k)J \u)l(p,k) wi(p,/c) ( 
I 23 (p, fc) = - 2e a ^^p^—- [ 1 




, ^2 (p, fe) w£ (p, fc) / \ w| (p, fc) wf (p, fc) 
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(127a) 



77 + -7777: . ( 127b ) 



^ 77 ^ 77777 I ( ~ 7777 7777 I » W 



^ + l —A . (127d) 



V ^3 (p, fe) (p, fc) J \ u\ (p, fc) u>\ (p, fc) 

We have taken into account that the still to be carried out loop integrals over d 3 p yield the 
same contribution for I 22 (p, /c,i) and 7 33 (p, A;,i), for Ii2(p,k,t) and Iis(p,k,t) and for l2i(p,k,t) 
and Isi(p,k,t). Hence, in each case, these contributions can be taken together to one single 
contribution, i.e., 

7 22 (p, k, t) = I 22 (p, k, t) + I 33 (p, k, t) = 2/ 22 (p, k, t) , (128a) 
/i 2 (p, k, t) = / 12 (p, k, t) + / 13 (p, k, t) = 2/ 12 (p, k, t) , (128b) 
hM k, t) = I 2i {p, k, t) + I 34 (p, k, t) = 2/ 24 (p, k, t) . (128c) 

The =-sign denotes that the equalities hold with respect to the integration over d 3 p. As the next 
step, we investigate the asymptotic behavior of the different Iij(p,k) for p — > 00 to determine 
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whether the integration over the loop momentum is finite, 



I U (P, k) =e 2 (l + x 2 ) (m ' ^f c? + O ( 1 ) , 
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r r^M e 2 (m c - m b ) 3 (l - x 2 )(l - x) , n fl\ noQ , 

hiip, k) =— rr — — ^- + O [-?] . (129g) 



2coj:p 3 [p 2 (l - x 2 ) + m 2 ] \fi 

As the integration measure, d 3 p, still contributes another factor of p 2 to the integrand, the latter has 
to be of the order of 1/p 4 for the loop integral to be finite. We see, however, that the contributions 
describing quark (i = 2) or antiquark (i = 3) bremsstrahlung and the interference between these 
two processes feature terms decaying as 1/p 2 and 1/p 3 in each case and that these terms do not 
cancel each other. Thus, the overall integrand behaves as 1 jp 2 for large p, which means that the 
loop integral is linearly divergent. 

In order to handle this divergence, we note that l22(??) k) and l23(p, k), from which this divergence 
arises, scale with the Bogolyubov particle number oc 1/p 2 for large p. This behavior is an artifact 
from the instantaneous mass shift [3^, [3^| (cf. section IIIip and can be regulated if the mass shift 
is assumed to take place over a finite time interval, r. This will be confirmed below in section 
IIVBI Hence, from the conceptual point of view, the linear divergence in the loop integral does not 
require a renormalization and is here regulated by cutting the loop integral at p = Ac- 



2. Asymptotic photon spectra 

Fig. [TBI shows the resulting photon spectra for different values of Ac. As in 0, [3^, we have 
chosen m c = 0.35 GeV and m b = 0.01 GeV. One can see that the photon spectrum drops as 
l/a;| in the ultraviolet domain such that the total number density and the total energy density of 
the emitted photons are logarithmically and linearly divergent, respectively, for given Ac- We will 
investigate below if this is also an artifact of the instantaneous mass shift. 

Furthermore, we investigate the dependence of the photon yield on the values of the constituent 
mass, m c , and the bare mass, m b , which is depicted in Fig. [TJl For small mass shifts, the photon 
yield scales with (m c — m b ) 2 , which is particularly visible for large Ac. This can be understood by 
taking into account that the dominant contributions to (|123p are given by (j!26b[) and (|127cf) , which 
for large p both scale with occupation numbers oc (ra c — m b ) 2 jp 2 . For large mass shifts, i.e., for 
m b —> at fixed m c , the photon yield starts to deviate from this quadratic scaling. In that limit, 
the photon yield diverges due to a collinear and an anti-collinear singularity in the loop integral. 
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FIG. 13: Photon spectra for an instantaneous mass shift for different values of Ac- In each case, they decay 
as l/uji in the ultraviolet domain. 
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FIG. 14: Dependence of the photon yield on the magnitude of the instantaneous mass shift for cvr = 1.0 GeV 
(left panel) and = 5.0 GeV (right panel). For small mass shifts, the yield scales oc (m c — mj) 2 , which is 
particularly visible for large values of Ac. 



This issue will also be discussed in greater detail in appendix IF1 From the phenomenological point 
of view, however, it is not a serious problem since the quark masses stay finite even in the chirally 
restored phase. 



3. Summary of results 

So far our investigations on chiral photon production have shown that the scenario of an in- 
stantaneous mass shift essentially comes along with three unphysical artifacts. 

Firstly, the loop integral entering (|62p features a linear divergence. In particular, this divergence 
arises from the contributions describing quark and antiquark bremsstrahlung and the interference 
between these two processes. It is caused by the quark and antiquark occupation numbers scaling 
oc (m c — rrib) 2 /2p 2 for p 3> m;,,m c . The latter is an artifact of the instantaneous change of the 
quark mass, and the mentioned divergence has been regulated by a cutoff at p = Ac- 

Furthermore, we have seen that then the asymptotic photon spectra decay as l/u;2 for any 
fixed value of Ac, which means that the total number density and the total energy density of the 
produced photons are logarithmically and linearly divergent, respectively. 

Finally, the asymptotic photon yield diverges for mj — > 0. We have demonstrated that this 
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divergence is due to a collinear and an anticollinear singularity in the loop integral over the con- 
tributions (|126bp and (|126cp describing quark and antiquark bremsstrahlung and quark-antiquark 
annihilation into a photon, respectively. This is, however, a less serious problem than the two 
previous ones as it can be circumvented by leaving the bare mass, m&, finite. The latter is justified 
from the phenomenological point of view since the quarks masses stay finite even in the chirally 
restored phase. 

Hence, as the next step, we have to determine if and to which extent these problems are regulated 
if the chiral mass shift is assumed to take place over a finite time interval, r, which corresponds to 
a physically more realistic scenario. 



B. Mass shift over a finite time interval 



1. Calculation of photon numbers 



For the general case of a mass shift over a finite time interval, r, both the time evolution of 
the fermionic wavefunctions, ipp 8 ,-\i( x ) j an d the time integrals entering (|59p require a numerical 
treatment. Hence, we have to find a way to extract the physical photon numbers from (|59p . i.e., the 
contributions which persist after taking the successive limits t — > oo and e — > 0. For this purpose, 
we consider the photon self-energy in terms of positive- and negative-energy wavefunctions again, 



in£,(Mi,f 2 ) = <?Y. 



d 3 p 



(2tt) 3 



(130) 



j f t{p,k,t 1 )j*(p,k,t 2 ) , 



(2tt)3- 

where we have introduced the effective current 



3ii(p, k, r, s, t) = ei/i'-^t)^ 



(t) 



(131) 



Next, we rewrite 



= 3™ ST (P, k, r, s, t) + jl(p, k, r, s, t) , (132) 
In analogy to V4 r -H (*)> t ne expressions ip^ (t) are defined according to the relation 



P,r,U 

Moreover, we have introduced 

•MST/-. 



(133a) 
(133b) 



in the second step. Eq. (|133ap vanishes if the quark mass is kept at its initial constituent value, 
m c , for all times, t, but (|133b[) does not. As a consequence, these expressions can be considered as 
mass-shift (MST) and vacuum contribution to ()131|) . With the help of (|133p . we can decompose 
the photon self-energy according to 



in< (k, h, t 2 ) = m<>°(k, t u t 2 ) + in^ MST (^, t u t 2 ) + in<^ INT (fc, t u t 2 ) 



(134) 
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with the individual contributions given by 

d 3 p 



r,s 

r,s ^ 

in<^ NT (fc,ti,t 2 ) = E 



f u (p, k, r, s, h)$*(p, ft, r, s,t 2 ) , 
d 6 p , 



(2tt) 3J 

3 



j^(p,k,r,s,h)j^(p,k,r,s,t 2 ) 



(27T) 

d 3 p 
(2^)3 



(135a) 
(135b) 



j"f ST (p,k,r,s,h)£'*(p,k,r,s,t 2 ) 



+j°Jp,k,r, S ,t 1 )j™ ST '*(p,k,r, S ,t 2 ) 



(135c) 



Expressions ()135ap and (|135bp describe the contributions from the vacuum polarization (VAC) and 
the mass shift (MST) to the overall photon self-energy, respectively, whereas (|135cp characterizes 
their interference among each other (INT) . Hence, it is convenient to split up (|59[) accordingly 



J2« - • 



k d 3 xd 3 A; 



d e n e Jt) 



k d 3 xd 3 A; 



d 6 n^(t) 



VAC 



+ 2w fc , 3 , 3; 



+ 2w 



d 6 n^(i) 



MST 



k d 3 xd 3 k 



(136) 



INT 



The individual contributions read 



d & nUt) 



2uj 



k d 3 xd?k 



2u) 



d 6 n^t) 



2uj 



k d 3 xd 3 ft 
d 6 n'(*) 



k d 3 xd 3 k 



VAC 



MST 



INT 



(2tt) 3 
7^ (ft) /•* 



dh I dt 2 m<£(k,t 1 ,t 2 )e 



<firt + . +_\„i^(ti-t 2 ) 



(2vr) 3 



dh 



dt 2 in<' MST (Mi,i2)e 1£ ^ 



W(tl-t 2 ) 



dh I dt 2 iU<^ mT (k,h,t 2 )^ { 



(137a) 
(137b) 
(137c) 



In appendix [A] it is shown that the contribution from the vacuum polarization (|135ap vanishes 
under the successive limits t — > oo and then e — > 0. We now demonstrate that the contribution 
from the interference term (|135cp is also eliminated by this procedure so that only the mass-shift 
contribution (|137bp remains and thus describes the physical photon number. For this purpose, we 
first rewrite the asymptotic contributions from (j!35bj) and (|135c|) at finite e by interchanging the 
time integrals with the loop integral over d 3 p. This leads to 



d b nl 

2uj ■ 



2uj 



k d 3 xd 3 k 
d 6 nf 



MST 



(2vr) 3 ^ 
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k d 3 xd 3 k 



j^(k) 



(2vr 

INT ' r,s 



d 3 p 
(2^)3 

d 3 p 
(2tt) 3 



I £ /(p,k,r,s)If,{p,k,r,s) 



Iu*(P,k,r,s)JC(p,k,r,s) 



(138a) 



+ Ja*(p,k,r,s)I*(p,k,r,s) 



In order to keep the notation short, we have introduced 



Jft&k, r,s) 



dtMt)ju ST (p,k,r, S ,t)e^ , 



-oo 
oo 



dtfS)fM k ^ s ^y u%t 



(138b) 

(139a) 
(139b) 
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The time integral in (|139b|) evaluates to 

J»(P, k, r, s) = eu c (p, s)^^v c {p + k, s)— ^ r , (140) 

e z + u{ (p, k) 

with uj\(p, k) given by (|118ap . To handle the time integral entering (|139ap . we first split 

/T poo 
dtf £ (t)j™ ST (p, k, r, s, t)e^ + / dtf £ (t)j™ ST (p, k, r, s, t)^ , (141) 
-oo JT 

where T 3> r. Next we take into account that for t > T, the fermionic wavefunctions have 
essentially turned into superpositions of positive- and negative-energy states with respect to the 
final bare mass, m b , i.e., 

4,*,t(*) = «^Ct(*) + Pi$i,.Jt) > ( 142a ) 

^.dt) = ^is^+^i,^) . ( i42b ) 

with the coefficients not depending on time. We have introduced the ~ notation in order to highlight 
that the expansion coefficients are generally different from those for an instantaneous mass shift 
given by (|95ap . With the help of (|142p . expression (|14ip is further evaluated to 

i;(p,k,r,s)= f tof&^f^ipXr^t)^* 

J —oo 

( _ e -[e-iw}tf$]T 

+e < a*fj f+j :u b (p, r) ltl v b {p + k, s) , 



+ a *p$p+k u b(Pi rh^u b (p + k,s) 



£ — l(jJ\(p, k) 
e -[ £ +io; 2 b (p,fc)]r 



+ P*plp + k v b{P, rh^v b {p + k, s) 



e + ioo%(p, k) 

c 



e— iojg (p,k)\T 



(143) 



e — iu]%(p, k) 



e" |e 



+ P*p 5 s+% v b^ r)^u b (p + k,s)- 

_ e -H<4(^)] T> | 

- « C (P, r)j^v c (p + fe, s) — tt- > • 

£-l<4(p,fc) I 
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Since the frequencies ()118a|) - (|118d[) are either positive or negative definite, taking the limit e — > 
leads to 

f-T 



I* (p, k, r, s) -> J„ (p ; k,r,s) = [ dt j^ ST (p, fc, r, s, t)e^ 



+ie 



®*jf(p + k u b{P, r )l^ h {p + ft*, s 



c 



-i^(p,fc) 



e iw|(p,fe) 



(144) 



w|(p, A;) 
i^(p,fe) 



g —4 



w 4 (p, ft) 



-» e 
u c {p,r)'~i^v c (p + k,s" 



We are allowed to interchange the limiting process with the remaining time integral since 
j^ ST (p, k, r, s, t)e lw fc* vanishes for t — > — oo and is thus integrable by itself on the time interval 
(— oo;T]. As we also have from (|140p 

J£(p, k, r, s) = eu c {p, r^v^p + k, s) - ^ + C(e 2 ) , for e -»• , (145) 

u>{ (p, ft) 

the interference contribution vanishes in that limit. Furthermore, the mass-shift contribution, 
which as a consequence of the above describes the actual photon number, turns into 

2 ^¥^¥k = WiF ^ / Wf 1 *^ ^ r ' s)Iv{ ^ ^ r ' s) ■ (146) 



It follows from (|146p and ()144p that solving the equations of motion (j36j) numerically on the time 
interval [— T; T] essentially provides all the information required to evaluate I^(p,k,r,s) and the 
asymptotic photon numbers (|146j) . We have j^f ST (p,k,r, s,t) fa for t < —T since T 3> r. We 
thus can approximate 

dt jjf ST (p, ft*, r, s, t)e^ « / T dt i^ ST (p, ft*, r, a , t)e^' . 

Hence, the numerical solution of (|36p on [— T; T] allows us to evaluate the time integral entering 
(|146p with sufficiently high accuracy. Based on this solution, the asymptotic expansion coefficients 
can be projected out from (|142p at t = T. 

So far, we have restricted ourselves to the scenario where the quark mass is only changed from 
its constituent value, m c , to its bare value, m&. When considering the second scenario, where the 
quark mass is first changed from m c to and then back to m c , the asymptotic photon number is, 
however, determined mostly in the same way. The only difference is that T then has to be chosen 
such that T S> r + with tl denoting the lifetime of the chirally restored phase. Moreover, 
one has to take into account that the fermionic wavefunctions have turned into superpositions of 
positive- and negative-energy states of mass m c instead of mass for t > T, i.e., 

4,*,t(*) = a p^ii,t(*) + wisA® ' (i47a) 

= WLft) + W C L^) ■ (i47b) 
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Accordingly, expression (|144p is replaced by 



T 



Ijp, k, r,s)= dt j^ bL (p, k, r, s, t)e^ 



+ie 



a 



iTifl+k ~ 1 ) r^vdp + k, s) 



e iwJ(p,fc) 

uf{p,k) 



- apSg+j&cip, r)^u c {p + k, s) 

e M(p^) 

+ (3tf% + kVc(p, rh^v c (p + fc, s) 



p 'p+fc 

-Pjfif+jpcip, r)lfj,u c (p + fe, s) 



u%(p,k) 

e -i^|(p,fc) 

Lol(p,k) 



For completeness, we mention that with the help of 



(fc) = J>^*(MK(M) 



(148) 



expression (|146p can be brought into the following alternative absolute-square representation 



2co 



d 6 n 



: d s xd s k (2tt) 



7r^3 



A,r,s 



d 3 p 



(27T) 



(149) 



Thus, the photon number is positive (semi-) definite and cannot acquire unphysical negative values. 
Furthermore, it vanishes if no mass shift takes place at all since we then have j^ ST {p, k, r, s, t) = 0. 



2. Numerical investigations and results 

First of all, we have to determine whether the linear divergence in the loop integral entering 
expression (j62[) for the photon yield is cured if the mass shift is assumed to take place over a 
finite time interval, r. For this purpose, we consider the cutoff dependence of the asymptotic 
photon number for different photon energies, ur, and different mass parameterizations, rrii(t), 
which is depicted in Fig. [15j As mass parameters, we have again chosen m c = 0.35 GeV and 
m b = 0.01 GeV. 

We see that the linear divergence, which has shown up in the loop integral for an instantaneous 
mass shift, is absent for both parameterizations rri2(t) and m^{t). In particular, the order of dif- 
ferentiability of the considered mass parametrization, rrii(t), is crucial for the saturation behavior 
of the loop integral. For m,2(t) being continuously differentiable once, the loop integral saturates 
at Ac — 10 GeV whereas it exhibits a much faster saturation already at Ac ~ 2.0 — 3.0 GeV for 
ms(t) being continuously differentiable infinitely many times. Since the latter parametrization de- 
scribes the most physical scenario, chiral photon production can be considered as a low-momentum 
phenomenon. 

As the loop integral is finite for a mass shift over a finite time interval, r, we can now turn to the 
UV behavior of the resulting photon spectra. Fig. [16] compares the resulting photon spectra for the 
different mass parameterizations. For m,2{t) and m^{t), a transition time of r = 1.0 fm/c has been 
assumed. Analogously to the particle spectra investigated in section [TTTl we see that the asymptotic 
photon spectra exhibit a strong sensitivity to the order of differentiability, i.e., the 'smoothness' 
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FIG. 15: Cutoff dependence of asymptotic photon numbers for 7712(f) (left panel) and 7713(f) (right panel). 
In both cases, we have chosen r = 1.0 fm/c. The parameterizations rrii(t) are shown in Fig. [5] and defined 
in The loop integral is rendered finite and its saturation behavior crucially depends on the order of 

differentiability of the considered mass parametrization. 
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FIG. 16: Asymptotic photon spectra for the different mass parameterizations given by Eqs. If one 

turns from an instantaneous mass shift (mi(t)) to a mass shift over a finite time interval ^^(t)), the 
photon spectra are rendered integrable in the ultraviolet domain. Furthermore, their decay behavior there 
is highly sensitive to the order of differentiability of the considered mass parametrization, rrii(t). 



of the considered mass parametrization, rrii(t). In particular, the decay behavior in the ultraviolet 
domain is suppressed from oc 1 /uj'i to oc 1 /lj® if we turn from mi (f ) (discontinuous parametrization) 
to 7772(f) (parametrization being continuously differentiable once). The logarithmic and linear UV 
divergences in the total photon-number density and the total energy density, respectively, are thus 
cured. Furthermore, if we consider the photon spectra for 7773(f), which is continuously differentiable 
infinitely many times and hence describes the most physical scenario, the photon numbers in the 
UV domain are suppressed even further to an exponential decay. 

As one can infer from Fig. [T71 the decay behavior of the photon spectra is highly sensitive to 
the considered transition time, r, for both 7712(f) and 7773(f). I n each case, the suppression of the 
photon numbers compared to the instantaneous case is the stronger the more slowly the mass shift 
is assumed to take place. As expected, both parameterizations reproduce the photon spectra for 
an instantaneous mass shift in the limit r — > 0. 

It should be mentioned that because of finite machine precision, it was first difficult to resolve 
the photon numbers numerically for 7773(f) an d r = 1.0 fm/c in the domain uj^ > 2.5 GeV. For that 
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FIG. 17: Asymptotic photon spectra for different transition times, r, for 7712(f) (left panel) and 7773(f) (right 
panel). The suppression of the photon numbers with respect to the instantaneous case is the stronger the 
more slowly (r increasing) the mass shift is assumed to take place. 



reason, the photon numbers in that domain have been extrapolated from those for 1.5 < lu? < 2.5 
GeV by performing a linear regression of the logarithms of the photon numbers. 

We shall briefly point out why we have to go to comparatively small transition times of r ~ 
0.02 fm/c for the photon numbers to be again of the same order of magnitude as for an instantaneous 
mass shift. The main reason is the different convergence behavior of the loop integral for different 
mass parameterizations. For the case of an instantaneous mass shift it features a linear divergence, 
which means that all momentum modes with p S> m c , nib an d p S> u)r contribute more or less 
equally to (|62p . For the case of a mass shift over a finite time interval, r, however, the loop integral 
is UV finite so that the contributions from the different momentum modes are suppressed with 
increasing p. This implies that the suppression of (|62p for given r with respect to the instantaneous 
case, m\(t), is the stronger the larger the value of Ac is chosen. Therefore, the larger Ac is chosen 
the smaller r has to be taken in order to approximately reproduce the photon yield for m\(t). 
This can also be inferred from Fig. [18] displaying the cutoff dependence of the photon yield for 
ur = 1.0 GeV and different transition times, r. 
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FIG. 18: Cutoff dependence of the photon yield for 7772(f) (left panel) and 7773(f) (right panel) at different 
values of r. For given finite r the suppression compared to the instantaneous case is the stronger the larger 
Ac is chosen. 



For completeness, we also investigate the dependence of the resulting photon spectra on the 
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magnitude of the mass shift, which is depicted in Fig. fl9l Similarly to the instantaneous case, the 
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FIG. 19: Dependence of the photon numbers on the considered bare mass, mj, for rri2(t) (left panel) and 
m^,(t) (right panel). As in the instantaneous case, the photon numbers increase with the magnitude of the 
mass shift and the behavior in the limit m& — > indicates a divergence therein. 

photon yield arising from the chiral mass shift increases with the magnitude of the latter. The 
change in curvature which appears for m& — > indicates a possible divergence in this limit which, 
in analogy to the instantaneous case, could arise from a collinear and/or anticollinear singularity 
in the loop integral entering (|62|) . This still requires further investigation. 

As for the asymptotic quark/antiquark occupation numbers, we also consider the scenario where 
the fermion mass is first changed from m c to toj and then back to m c to take into account the finite 
lifetime of the chirally restored phase. Fig. [20] shows the photon spectra for different values of Ac 
for both mass shifts taking place instantaneously which is described by fh\(t). We have assumed 
a lifetime of tl = 4.0 fm/c for the chirally restored phase. 
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to k [GeV] 

FIG. 20: Asymptotic photon spectra for both mass shifts taking place instantaneously. The loop integral 
in ()62l) is again linearly divergent. Furthermore, the photon spectra decay as if this divergence is 

regulated via a cutoff at p = Ac- 

As for the first scenario, the loop integral entering (|62|) exhibits a linear divergence. If this 
divergence is regulated via a cutoff at p = Ac, the resulting photon spectrum again decays oc 
in the ultraviolet domain and is hence not integrable. As one would expect from the first scenario, 
however, these pathologies are again artifacts from the (unphysical) instantaneous mass shifts and 
are resolved if both mass shifts are assumed to take place over a finite time interval, r. 
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FIG. 21: Cutoff dependence of asymptotic photon numbers for m,2(t) (left panel) and fh^{t) (right panel). 
In both cases, we have chosen r = 1.0 fm/c and tl = 4.0 fm/c. The parameterizations fhi(t) are shown in 
Fig. [5] and defined in (|105p . Like in the first scenario, the loop integral is rendered finite and its saturation 
behavior crucially depends on the order of differentiability of the considered mass parametrization. 



In particular, one can infer from Fig. [21] that the loop integral again saturates around Ac* = 
10 GeV and Ac = 2-3 GeV for rfi2(t) and rh^t), respectively. Moreover, Fig. I2"2l shows that 
the resulting photon spectra exhibit the same sensitivity to the order of differentiability of the 
considered mass parametrization, rrij(i), just like in the previous case. The photon spectra decay 
oc l/w§ for rfi2{t) (continuously differentiable once) in the ultraviolet domain and are suppressed 
further to an exponential decay for rh%(t) (continuously differentiable infinitely many times). 
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FIG. 22: Asymptotic photon spectra for the different mass parameterizations given by Eqs. (|105|) . As 
expected, the photon spectra are again rendered integrable in the ultraviolet domain if both mass changes 
take place over a finite time interval, r. Furthermore, they show the same sensitivity to the order of 
differentiability of the considered mass parametrization, rhi(t), as for a single mass change. 



As it must be, the photon spectra show the same sensitivity to the change duration, r, as for 
the first scenario, i.e., the suppression of the photon numbers compared to the instantaneous case 
is the stronger the more slowly that mass changes are assumed to take place. Moreover, both ifi2(t) 
and fh^{t) reproduce the photon spectra for the instantaneous case in the limit r — > 0. This is 
shown in Fig. [23J. 

Hence, we have seen so far that the general dependence of the photon numbers in the ultraviolet 
domain on the order of differentiability of rh(t) and the transition time, r, is the same as for the 
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FIG. 23: Asymptotic photon numbers for different transition times, r, for rfi2{t) (left panel) and fh^{t) (right 
panel). As in the previous scenario, the suppression of the photon numbers with respect to the instantaneous 
case is the stronger the more slowly (t increasing) the mass shifts are assumed to take place. 



first scenario, which one would also expect intuitively. Nevertheless, there are some differences in 
the dependence on the magnitude of the mass shift, m c — nib, which can be seen in Fig 1241 
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FIG. 24: Dependence of the photon numbers on the considered bare mass, mj, for 7712 (i) (left panel) and 
m3(/j) (right panel). Like for the scenario where the mass is solely changed form m c to m^, the photon 
numbers increase with the magnitude of the mass shift, m c — nib. Their behavior in the limit m& — > does, 
however, not indicate a divergence. 



As in the previous case, the photon yield increases with the magnitude of the mass difference for 
given photon energy, cj^. The crucial difference, however, is that the curvature does not change for 
TTife — > 0. This indicates that, in contrast to the first scenario, the photon numbers converge in that 
limit and that the loop integral entering (|62p does not feature a collinear and/or an anticollinear 
divergence therein. 

We have seen in section IIIII that the asymptotic quark and antiquark occupation numbers for 
p S> m c , mi, are modified by a factor of 4 sin 2 jjtx when turning from rrii(t) to rhiit) (i = 1, 2, 3) for 
any given transition time, r. In contrast, the asymptotic photon numbers in the ultraviolet domain 
do not exhibit a similar modification by a factor of 4 sin 2 co^tl ■ Solely for 777,3 (t) the photon spectra 
exhibit a slightly oscillating behavior (see e.g. right panel of Fig. 123]) in the photon momentum, k, 
for sufficiently large values of r. For 777,3(7;) and r = 1.0 fm/c this behavior can only be displayed 
up to = 2.5 GeV since the photon numbers for cj^ > 2.5 GeV have again been extrapolated 
from those for 1.5 <wr< 2.5 GeV by a linear regression. 
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Even though one might expect a similar modification as for the asymptotic quark/ antiquark 
occupation numbers in the first place, there are two important aspects to be taken into account. 
On the one hand, the asymptotic photon numbers incorporate the entire history of the fermionic 
wavefunction and hence of the quark and antiquark occupation numbers, which are extracted from 
the former. In particular, the occupation numbers partially coincide with the asymptotic ones 
of the first scenario between the two mass shifts. On the other hand, it follows from (j62|) that 
the dependence of the wavefunction parameters on the fermion momenta is integrated out when 
determining the asymptotic photon numbers. Upon this procedure, a possible oscillating behavior 
in the individual contributions to (|62p from the different momentum modes can get lost again. 

The latter aspect is supported when comparing the asymptotic photon numbers for both param- 
eterizations rrii(t) and rhi(t), which is done in Fig. [25] for fixed photon energies, oj^, and different 
magnitudes of the mass shift, m c — mf,. There the dotted lines represent the photon spectra of 
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FIG. 25: Comparison of photon numbers for m^ify and m^it) (left panel) and for m^t) and m-s{t) (right 
panel). For small mass shifts, m c — mj,, the photon numbers differ by a factor of roughly 2, which is 
particularly distinctive for mi(t) and rri2{t), respectively. 



rriiit) multiplied by a factor of 2 in each case. Hence, for small magnitudes of m c — mb, the photon 
numbers roughly double if the quark/antiquark mass is switched back to its constituent value, m c . 
The latter feature is particularly distinctive for 7712 (£)■ Such a result is understandable since both 
mass shifts are expected to give a comparable contribution to the asymptotic photon yield (|62|) . 
But in particular, integrating out an additional factor of 4 sin 2 pr^ gives rise to an overall rescaling 
by a factor of roughly 2 if the integrand does not change significantly over the periodicity interval 
Ap = tt/tl- For increasing m c — m&, the asymptotic photon numbers for m,i{t) and fhi(t) start to 
deviate from this ratio as the different scaling behavior for mj, — > starts to manifest itself. 

In contrast to 3CH32f| . our asymptotic photon numbers arising from first-order QED processes 
are UV finite for both scenarios if the mass shifts are assumed to take place over a finite time 
interval, r. Hence, it is convenient to compare them to leading-order thermal contributions. We 
note again that first-order QED contributions vanish in a static thermal equilibrium. There the 
first non-trivial contribution starts at two-loop order. Since a loop expansion does not coincide 
with a coupling-constant expansion, resummations are necessary to obtain the thermal rate at 
quartic order in the perturbative coupling constants, i.e. at linear order in a e and at linear order 
in a s 50(]. Fig. [26] shows the photon numbers from first-order chiral photon production together 



with leading-order thermal contributions. The latter have been obtained by integrating the leading 
order thermal rates taken from [5(3] over the lifetime interval of the chirally restored phase, tl = 4.0 
fm/c, at a temperature of T = 0.2 GeV. 
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FIG. 26: Comparison of first-order mass-shift contributions to time-integrated thermal rates for rrii(t) (left 
panel) and fhi(t) (right panel). For i = 3 most likely characterizing a physical scenario, the photon spectra 
arising from the chiral mass shift are clearly subdominant for lo^ > 1.0 GeV. 



In this context, the photon spectra for rhi(t) are the more meaningful ones, as the underlying 
scenario considers the finite lifetime of the chirally restored phase during a heavy-ion collision. 
Nevertheless, we see that for both m^(t) and m^(t), which are continuously differentiable infinitely 
many times and thus represent the most physical scenario, the photon numbers for ojj; > 1.0 GeV 
are subdominant compared to those obtained from (time-integrated) thermal rates for phenomeno- 
logically reasonable choices of r and T. 

To summarize, we have seen that, if we turn from an instantaneous mass shift to a mass shift 
over a finite time interval, r, the linear divergence in the loop integral is regulated. Furthermore, 
the asymptotic photon spectra are integrable in the ultraviolet domain if the time evolution of 
the quark masses is described in a physical way. Finally, the decay behavior shows a strong 
sensitivity to the considered transition time, r, of the quark mass and we recover our results for an 
instantaneous mass shift as r — > 0. The dependence on rh(t) and r is the same if the quark mass 
is also restored to its constituent value, m c , as to mimic the finite lifetime of the chirally restored 
phase. 

In particular, for mass parameterizations that are continuously differentiable infinitely many 
times and thus represent the most physical scenario, our photon numbers are subdominant with 
respect to those arising from integrated thermal rates in the UV domain for a physically sensible 
transition time, r, and temperature, T. 

Nevertheless, the dependence of the photon number on the bare mass, mj,, indicates that the 
loop integral entering expression (|62p for the asymptotic photon numbers still features a collinear 
and/or an anticollinear singularity in the limit mj —> for the first scenario (mass solely changed 
from m c to m&) whereas there is no such indication in the second one (mass changed back to m c ). 
Even though this aspect still requires further investigations, it is important to point out again that 
a possible divergence in the massless limit can be circumvented leaving m& finite. 



C. Summary of results 



We have seen that our prescription of chiral photon production eliminates possible unphysical 
vacuum contributions and leads to photon spectra being integrable in the UV domain for physical 
mass-shift scenarios. The crucial difference compared to been the consideration of 

asymptotic photon numbers. For this purpose, we have switched the electromagnetic interaction 
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adiabatically according to (|44p and determined the photon numbers for t — > oo. Only at the end 
of our calculation, we have taken e — > 0. It shall be stressed again that adhering to the correct 
order of limits is indeed crucial for two main reasons. 

First, the interpretation of (|59p as a photon number is only justified in the limit t — > ±00 for 
finite e where the electromagnetic field is non-interacting. As a consequence, taking first e — > at 
some finite time, t, is questionable, as we then would have an interacting electromagnetic field such 
that the interpretation of (|59p as photon number is not justified. Moreover, such an interpretation 
remains doubtful for t — > 00. Since we would have taken e — > before, the electromagnetic field 
would not evolve into a free one for t — > 00. A similar problem occurs when only using an adiabatic 
switching-on of the interaction for t — > —00 but no adiabatic switching-off for t — > 00. Such a 



procedure has been suggested in [5J| in order to implement to implement initial correlations at 
some t = to evolving from an uncorrelated state at t — > —00. 

Second, interchanging both limits comes along with a violation of the Ward-Takahashi identities. 
To see this, we consider the unphysical scenario, where the electromagnetic interaction is switched 
on an off again instantaneously at t = =Fl/e, i.e. 

H EM (t) -> e £ (t)H EM (t) , with 9 e {t) = e(^-t 2 ^j and e > . (150) 

We shall show in greater detail in appendix [C] that if we replace f £ (t) by 9 £ {t) in (|59p and consider 
this expression for free asymptotic fields, i.e., if we first take t — > 00 and then s — > 0, we obtain the 
same result for (flTS7b"l) as we would also have obtained if we had adhered to (|44p but interchanged the 
limits instead. Consequently, interchanging the limits for (|44p is formally equivalent to switching 
the electromagnetic interaction on and off again instantaneously, which is unphysical. 

In order to emphasize this aspect, we consider again the Ward-Takahashi identities ()61|) for 
the photon self-energy. We have shown in section |TT] and also demonstrated in greater detail in 
appendix [B] that these identities are fulfilled when coupling the quark fields to a time-dependent 
scalar background field, (ft(t). But it follows from (|59p that a switching of the electromagnetic 
interaction leads to a modification of the one-loop photon self-energy, i.e., 

in^Mi.ta) m<f{k,h,t 2 ) = / e (ti)/ e (t 2 )in^(^ti,t2) • (151) 

The electromagnetic interaction is then fully persistent at any given time, t. Because of the 
additional factor of f e (ti), relation (|6ip is not fulfilled in the first place, but we instead have 

d tl m<£(k,h,h) - ik j mff(k,h,t 2 ) = -/ e (*i)/ e (* 2 )iiiSi B (fc,ti,t2) • (152) 

The time derivative of f e (ti), however, reads 

f £ (t) = e sign(i)e- £ l*l . (153) 

As a consequence, the r.h.s. of (|152p vanishes for e — > and, accordingly, (|6ip is restored again. 
To the contrary, if we replace f e (ti) by 9 £ (ti), we have no such restoration since 

«.W = «(i-*)i(i + «)-«(i+*)i(i-«) . (154) 

Hence, switching the electromagnetic interaction as in (|150p leads effectively to a violation of the 
Ward-Takahashi identities. Since ()150p leads to the same asymptotic photon numbers as (|44p with 
the interchanged order of limits, we conclude that the latter procedure also leads to an effective 
violation of (EU). 
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From this point of view, the Ward-Takahashi identities actually imply constraints for a physically 
reasonable definition of photon numbers. Their violation implies an inconsistency of such definitions 
with £7(1) gauge invariance which, in turn, usually leads to unphysical artifacts. This issue has 
also been pointed out in the context of the Thomas-Reiche-Kuhn sum rules [52], which can be 
considered as direct consequences of charge conservation and gauge invariance of QED. It has been 
shown that these sum rules impose severe constraints concerning the applicability of trans por t 
approaches on photon production from non-equilibrated hot hadronic matter presented in 5^-57]. 



V. SUMMARY, CONCLUSIONS AND OUTLOOK 

In this work, we have investigated photon emission during the chiral phase transition in the 
early stage of a heavy-ion collision. During this phase transition, the quark mass is changed from 
its constituent value, m c , to its bare value, mp , which leads to a spontaneous non-perturbative pair 



production of quarks and antiquarks [38|, [39J . This effectively contributes to the creation of the 
QGP, and we have investigated the photon emission arising from this creation process. 

As inBfl, the change of the quark mass has been modeled by a scalar background field in the 



QED Lagrangian. We have restricted our considerations to first-order QED processes. Such pro- 



cesses have also been investigated in [4fj] for fermions coupled to a time-dependent electromagnetic 
background field [4l|]. They are kinematically allowed since the background field acts as a source 
of additional energy. For a proper treatment of the non-perturbative nature of the pair-creation 
process the coupling to the background field must be resummed to all orders. In order to achieve 
this, we have constructed an interaction picture including the coupling to the source field. The 
photon yield has then been obtained by a standard perturbative QED calculation. In the course 
of our calculations, we have essentially pursued an in/out description, where the photon numbers 
have been extracted in the limit t — > oo for free asymptotic states. In order to obtain such states 
we have introduced an adiabatic switching of the electromagnetic interaction, i.e., 

H E M(t) -> f £ (t)H BM (t) , with f e (t) = e- £ l*l and e > . (155) 

The photon numbers then have been considered for t — > oo and we have taken e — > at the very 
end of our calculation. 

Before turning to our investigations on photon production, we have first provided a digression 
on pair production of quarks and antiquarks. There, we have extended the investigations in 38|, 3^] 



to the time evolution of the quark and antiquark occupation numbers. We have restricted ourselves 
to the pair production arising from the chiral mass shift only and not taken into account radiative 
corrections. As in 3^, 39], we have compared different mass parameterizations 



m c + mj m c — m\, 
mi(t) = sign(t) , (156a) 

m 2 (t) = - sign(t) ( 1 - e- 2 '^) , (156b) 



m 3 (t) = b — — - tanh — . (156c) 



2 2 \t 

Here m c and mj, denote the initial constituent mass and the final bare mass, respectively. m\(t) 
describes an instantaneous mass shift at t = whereas rri2(t) and m^{t) denote a mass shift over 
a finite time interval given by r. 

We have seen that for the case of an instantaneous mass shift, the quark and antiquark occu- 
pation numbers for t > scale oc (m c — nib) 2 /p 2 for p 3> m c ,mfo, with p denoting the (absolute 
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value of) the fermion momentum. This means that the total particle number density and the total 
energy density of the fermionic sector are linearly and quadratically divergent, respectively. 

If we turn from an instantaneous mass shift to a mass shift over a finite time interval, r, the 
mentioned problems are resolved in the asymptotic limit t — > oo. For rri2{t) which is continuously 
differentiable once, the occupation numbers for p ^> m^m^ and p » 1/r are suppressed to 
oc {m c — m,b) 2 /p 6 T 4 such that both the total particle number density and the total energy density 
are rendered finite. Moreover, if we turn from m,2{t) to 7773 (i) which is continuously differentiable 
infinitely many times, the occupation numbers are suppressed further to an exponential decay. 

At finite times, t, however, the occupation numbers decay oc 1/p 4 for large p. This means that 
the total particle number density is finite whereas the total energy density features a logarithmic 
divergence. In this context, it is nevertheless important to point out that only asymptotic particle 
numbers describe observable quantities and that the analogous expressions at finite times do, in 



general, not allow for a similar interpretation 43]. 



For that reason, we have concentrated our investigations on photon production to free asymp- 
totic states. We have again compared the different mass parameterizations given by (|156ap - (|156cp . 
As one would expect from our investigations on particle production, where the total particle num- 
ber density is linearly divergent for an instantaneous mass shift (|156ap . the asymptotic photon 
numbers feature very similar pathologies for such a scenario. In particular, the unphysical artifacts 
we have encountered are the following: 

• The loop integral entering our expression for the photon yield features a linear divergence. 
We have shown that this divergence does, indeed, arise, from the decay behavior of the 
quark/antiquark occupation numbers oc (m c — m^/p 2 for p 3> m c ,mt,- As this decay 
behavior and thus the mentioned divergence is an artifact of the instantaneous mass shift, 
we have regularized the latter by cutting the loop integral at p = Ac- 

• After this regularization procedure the asymptotic photon spectra decay oc for large 

photon energies, uj? = |fc| (k denotes the three-momentum of the photon), and are thus not 
UV integrable. As in the total photon numbers density and the total photon energy 



density are logarithmically and linearly divergent, respectively. 

We then turned from an instantaneous mass shift to a mass shift over a finite time interval, r. 
This renders the loop integral finite. Moreover, the fermion-momentum range from which the 
main contributions to the photon numbers arise is highly sensitive to the order of differentiability, 
i.e., 'the smoothness', of the considered mass parametrization, rrii(t). For 777,2 (i) (continuously 
differentiable once), we have contributions from quark and antiquark momenta up to p ~ 10 GeV 
whereas we have contributions only up to p ~ 2.0 — 3.0 GeV for 7773 (i) (continuously differentiable 
infinitely many times). Since the latter parametrization describes the physically most realistic 
scenario, chiral photon production can accordingly be considered as a low-energy phenomenon. 

Furthermore, the resulting photon spectra are rendered integrable in the ultraviolet domain if 
the mass shift is assumed to take place over a finite time interval, r. In particular, their decay 
behavior in the ultraviolet domain also crucially depends on the smoothness of the respective mj (i) . 
For 777-2 (t), they decay as l/^l in the ultraviolet domain whereas they are suppressed even further 
to an exponential decay for m^{t). For both parameterizations, the suppression of the photon 
numbers with respect to the instantaneous case is the stronger the more slowly the mass shift is 
assumed to take place. As expected, both 7772 (t) and 7773 (t) reproduce the photon spectra for an 
instantaneous mass shift in the limit r — > 0. 

The asymptotic photon spectra, considered as the only observable ones, hence do show a very 
similar sensitivity on the mass parametrization, 777(t), as the asymptotic particle spectra, as to be 
expected. In particular, the logarithmic divergence in the energy density of the fermionic sector at 
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intermediate times does not appear in form of a similar pathology in the photonic energy density. 
This could not be excluded a priori since the asymptotic photon numbers incorporate the entire 
history of the fermion-wavefunctions and hence of the corresponding quark/antiquark occupation 
numbers. 

In order to take into account the finite lifetime of the chirally restored phase during a heavy- 
ion collision, we have also considered the scenario where the quark mass is changed back to its 
constituent value, m c , after some time interval, tl- Here we have again compared different mass 
parameterizations 



rhi(t) 
m 2 {t) 



m c + nib 



m c + nib 



sign \-f-f 

—sign [t+Y 
xsign ( — -t 



\2t+T L \/l 



\2t-T L \/l 



m c + m b m c -m b ^ ( 2t + r L \ { t l - 2t 

m^(t) = = tann | | tann 



(157a) 

(157b) 
(157c) 



The general dependence of the photon numbers on the 'smoothness' of the considered mass 
parametrization, rhi(t), and the transition time, r, has been the same as for the previous case. 
In particular, for mass parameterizations being continuously differentiable infinitely many times, 
our photon numbers have been subdominant compared to time-integrated thermal rates in the 
ultraviolet domain. 

A principal difference between both scenarios, with and without restoration of the mass to its 
initial value, has only been observed in the dependence on the magnitude of the mass shift, m c — nib. 
In the first scenario, this dependence indicates a divergence of the photon numbers in the limit 
nib —> whereas this is not the case for the second scenario. 

In conclusion, we have provided an ansatz for chiral photon production that eliminates possible 
unphysical contributions from the vacuum polarization and, furthermore, renders the resulting 
photon spectra UV integrable if the time evolution of the quark mass is modeled in a physically 
realistic manner. In contrast to 1371 ] . the photon numbers have not been considered at finite 

times, t, but for free asymptotic states. In this context, we have seen that a consistent definition 
of photon numbers is actually only possible for free asymptotic states. A similar interpretation 
of the respective expression is usually not justified for finite times, t, as we then do not have free 
asymptotic states. The same problem occurs if the electromagnetic interaction is only switched on 
adiabatically from t — > — oo but not off again for t — > +oo. This procedure has been suggested in 
511 ] such as to implement initial correlations at some t = to evolving from an uncorrelated state at 
► — oo. 

Moreover, the consideration of 'photon numbers' at finite times again comes along with a viola- 
tion of the Ward-Takahashi identities. Consequently, our investigations support the corresponding 
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concern raised in 33|, 34 ] towards 3(3-32]]. Our outlook to future investigations is hence as follows. 

Primarily, our results indicate that the problems with the divergent contribution of the vacuum 
polarization and/or the decay behavior of the resulting photon spectra appearing in 
result from an inadequate definition of the photon numbers for finite times, t. This, in turn, gives 
rise to the question whether they can be cured if the photon numbers are considered in the realm 
of asymptotic fields. In this context, the actual role of the Ward-Takahashi identities, which are 
violated in l37j |. but conserved by our approach on chiral photon production, still requires 

a more profound consideration. 

Furthermore, it is of particular interest whether our asymptotic definition of photon numbers 
can be extended to finite times in a consistent manner and/or which alternative quantities 
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can be considered to provide a proper real-time description of the electromagnetic sector during a 
heavy-ion collision. One promising candidate could be the gauge-invariant field strength tensor, 



and the quantities derived from it. One example is the energy density of the electromagnetic field, 
which is given by 



i(x) = F^{x)Fj{x) - - A ^F^{x)F Pa {x) 
=F^(x)F*(x) -\p^(x)Fp a {x) 



E l + B' 



fi=u=0 



with T]^ = diagjl, —1,-1,-1} and E and B denoting the electric and magnetic field operators, 
respectively. Investigations on such quantities are, of course, not restricted to chiral photon pro- 
duction but can also be based on the approach by Boyanovski et al. To obtain physically 
reasonable results in this approach might, however, require the consideration of correlated initial 
In 3CH32I]. the authors assume that the system is in thermal equilibrium with 



states 



5J, [5 



respect to the strong interactions, i.e., the initial state is specified by ([T]). The authors effectively 
neglected initial correlations characterized by the interaction part of Hqqd by using the one-loop 
approximation for the photon self-energy. It has, however, been shown that the evaluation of the 
photon self-energy cannot be performed by simple power counting in the strong coupling constant, 



a s , but instead requires a resummation of the so-called ladder diagrams [5CJ, I6CH63I] . Within this 



procedure, initial correlations are included, which leads to the question whether this regulates 
possible problems with the vacuum polarization and the UV behavior. 

As the problem of photon production as well as of the time evolution of the electromagnetic 
energy density can be reduced to the calculation of the photon self-energy, i.e., the current-current 
correlator, at first order in a e , another alternative is to address the question of finite-lifetime 
effects within the 2PI-approximation of the effective action. This has already been suggested in 
(35)]. Even though the conservation of the Ward-Takahashi identities still turns out to be difficult, 
the 2PI approach has the crucial advantage that it does not make any ad hoc assumptions about 
the two-time dependence of correlation functions. 
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Appendix A: Representation of photon yield as an absolute square 

In this section, we show that (|59p can be written as an absolute square of a first-order QED 
transition amplitude and is thus positive (semi-) definite. We also emphasize that keeping the 
correct sequence of limits leading to (|62|) is crucial to eliminate possible contributions from the 
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vacuum contribution (|64p to the photon self-energy and that an interchange comes along with a 
divergent contribution from it. In order to show that the photon yield can be written as an absolute 
square, we first undo the contraction 



m<(k,t 1 ,t 2 ) = ^(k)m<(k,t u t 2 ) 



Furthermore, we have 



in^(Mi,*2) = y {%q\jl^ k ^)ivA k M)\%q) 



^E<°^I$,j(M2)|/> </bV(Mi)|o 



ml ' 



(Al) 



(A2) 



with j^j{k,t) given by 



Jm,j(M) = d x j M ,j(x)e" 



-ik-x 



(A3) 



and |/) denoting an orthonormal basis of the fermionic Hilbert subspace. Together with (|48[) and 
i2j) . we can write (|59p as 



2u 



d 6 n^(t) _ 
fe 'd 3 xd 3 /t ~ (2vr) 3 F 



— y 

■\3T/ 



A,/ 



d«/ e ( U ) (/| 6?{%, A)i M ,j(fe, «) |0 OT -> e^ u 



(2vr) 3 F ^ 
v ; A,/ 



dS</|^(fc,A)^,j(x)|0 M -)e 



with the underline denoting 

[ d i x= f d 3 x f duf £ (u) with / e (t) = exp(-e|t|) . 

J J J —CO 

Using the relations (|22aft and (|25ap . expression ()A4|) can be rewritten further as 



d fa n^) _ 
d 3 xd 3 fc ~ (2it) 3 V 



— y 



A./ 



dW/;M 



= ^y 

y ' A,/ 

To keep the notation a little bit more shorthand, we have introduced 



-Hj(x) 



duf E (u)(f;k,X Hj(u) 



f;k,x) = \f) 



k, A 



|0) = |(W}®|0 7 ) . 



(A4) 



(A5) 



(A6a) 
(A6b) 



Hence, ()A5[> shows that (f59|) can be written as the absolute square of the (space-time integrated) 
transition amplitude at first order in e between the initial vacuum state and a final state containing 
a single photon of momentum, k, and polarization, A, which is summed over together with the 
fermionic degrees of freedom, |/). Furthermore, we point out that when taking successively the 
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limits t — > oo and then e — > 0, expression (|59[) does not contain any unphysical contributions from 
the vacuum polarization, i.e., Eq. (|59p vanishes if the quark mass stays at its constituent value, 
m c , for all times, t. In this case, we have 



in^Mi.taJe^* 1 -* 3 ) =iU< G {kM -i 2 )e 



\uAti-t-l) 



2e^ 



Accordingly, (|59p turns into 



d 6 n e Jt) 



f d3p Ji + pzQra + + 

7 (2tt)3 1 + E%E°_ , 



k d 3 xd 3 k (2vr) 3 
Letting t — > oo first, we obtain 

/oo 
-oc 

and hence 



P p+k 
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(A7) 
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Taking now the second limit £ — > 0, we finally get 
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££££ - 
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2s 



p+fc p 



=0 , 



(A9) 



(A10) 



where we have taken into account that EZ + ^ + E^+co^ is positive definite in the second step. Thus, 
we have shown that the vacuum polarization does not contribute to the asymptotic photon yield. 

For completeness, we also emphasize that keeping the correct order of limits, i.e., taking first 
t — >■ oo and then e — > 0, is indeed crucial to eliminate the contribution from the vacuum polarization. 
To clarify this, we show that taking e — > first at some finite time, t, and t — > oo afterwards comes 
along with unphysical artifacts. Such a procedure corresponds to switching on the electromagnetic 
interaction solely from — oo to some time, to, i.e., with the regulator function 

g e (t) = 0(to - *)e £ (*-*°) + 9(t - t ) , 



which has been suggested in [5 If ] to implement initial correlations at t = to evolving from an 
uncorrelated initial state at t — >■ — oo. In this case, the photon yield can still be transformed to 
an absolute-square of the form (|A5j) with f £ (t) replaced by g e (t). If one first takes e — > at some 
finite time, t, however, one encounters a divergent contribution from the vacuum polarization. In 
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order to see this, we again split the photon self-energy into the stationary vacuum polarization and 
a non-stationary-mass shift contribution, i.e., 



illy (A;, h,t 2 ) = in^ (k, h-t 2 ) + iAn^ (k, h,t 2 ) 



(All) 



where iAIIy (k, t\, t 2 ) vanishes for t±,t 2 < to- Accordingly, we decompose (foUj) as 



2u 
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k j3 



d*xd 6 k 
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+ 2uj 
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k d 3 axl 3 A; 



(A12) 



CMS 



with the contributions from the vacuum polarization (VAC) and from chiral mass shift effects 
(CMS) reading 
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dtx / dt 2 g £ (ti)g £ (t 2 )UI^ (k,ti - t 2 )e 



iwg(ii-<2) 



dh / dt 2 5£ (ti)g e (t 2 )iAII<(^t 1 ,t 2 )e 
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After interchanging both time integrations with the loop integral the contribution from the vacuum 
polarization at time t is evaluated to 
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with the time integral given by 



du g £ {u)e v p J = -. 

oo e + i [W -, + E%+uj? 
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Taking now the limit e — > at this point, the vacuum contribution does not vanish but instead 
turns into 
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(A16) 



Since the integration measure d 3 p adds another factor of p 2 to the integrand, the loop integral 
entering (|A16P is linearly divergent. As (|A16|) is furthermore constant in time, it does not vanish 
in the limit t — > oo. Hence, interchanging the limits t — > oo and e — > comes along with a divergent 
contribution from the vacuum polarization. 

Since (|A16p is constant in time and hence already fully present before any mass-shift contri- 
butions characterized by iAH^,(k,ti,t 2 ) can show up, one might still argue that (lATBj) can be 
uniquely identified with the virtual cloud of the vacuum and, as a consequence, be subtracted 
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(renormalized) from the overall photon yield (|59p . This would only leave the contribution (|A13bj) . 
The main problem with this argument, however, is that this contribution on its own is not an 
absolute square, which can be seen by rewriting it as follows: 
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to 

+2Re 



dti / dt 2 g e {t 1 )gS2)i^l<{k,t l My^ {tl ~ t2) 
J to 

I dti f° dt 2 < 7£ (t 1 ) 5e (i 2 )iAn<(fc,t 1 ,t 2 )e i ^(' 1 - i2 )j) . (A17) 

J to J ~oo J J 



Above, we haven taken into account that iAIIy (k, t±, t 2 ) = for ti,t 2 < to- We have, however, 
iAIIy (fc, ti, t 2 ) / if either ti > to and t 2 < to or *i < *o an d *2 > to- This gives rise to the second 
term in ()A17j) . For this term, we have also made use of 

iAn<(£,ti,t 2 ) = iAn<(M 2 ,ti; 

Since the overall contribution from iIT^ (&:,ti — t 2 ) from the three domains 

• t%, t 2 > t0 J 

• t\ > to and t 2 < to , 

• ti < to and t 2 > to , 

moreover vanishes in the limit e — > 0, we can rewrite ()A17p as 
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(A18) 



In contrast to (|59p . there is no contribution to (|A18[) from the domain where t\, t 2 < to- Therefore, 
it is not possible to take (|A17|) into an absolute-square representation of the form of (|A5|) . but 
with f e (t) replaced by g e (t). Hence, it is not guaranteed that the (renormalized) photon yield is 
positive (semi-) definite, and it can, in principle, acquire unphysical negative values. 



Appendix B: Ward-Takahashi identities for ill< (k,t, u) 



In this section, we will provide an explicit verification of the Ward-Takahashi identities for the 
photon self-energy, Ul^ u (k, t, u), which read 



d t iU<{k,t,u)-ik l iU<(k,t,u) = 0, 

d t m<(k,t,u) -ikniif^k^u) = o. 



(Bla) 
(Bib) 



For this purpose, we first express iH^ u (k, t, u) in terms of wavefunction parameters. It follows from 
(1571) that 



in wl/ (fc,t,«) 
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r,s 

E 
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(B2) 
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With help of (|35p . the matrix elements are evaluated to 

-/?i t (t)a j) . + ^ 4 (t) Xr (<? • %) eras , (B3b) 

where ep and denote the unit vectors in the directions of p and p+k respectively. Furthermore, 
we have taken into account that 



7i = -7 



-a 1 
a* 



for spacelike indices, i. The representation of iH^ u (k,t,u) in terms of wavefunction parameters is 
obtained upon insertion of (|B3a}) - (|B3bj) into (|B2p for the respective indices and carrying out the 
summation over the two spin indices, r and s. For both indices being timelike, i.e., fj, = v = 0, we 
obtain 

-* f d^p ( 

m£ (k,t,u) = 2e 2 J ^3 |o ! t t (t)Q!^ + ^(t)Q!* + ^ 4 («)a^ 1 -(u) + 

ff-W+ftr (B4) 
P -^+ [^t(*)^(^(«)^t(«) + 

Here we have made use of the identities 

\Xr (<? ■ u) (a ■ v) Xs \ 2 = Tr {(«? • u) 2 (a ■ v) 2 } = 2u 2 v 2 , (B5a) 

r,s 

y~] Xs (d ■ u) (a ■ v) XsS rs = Tr{(<T • u) (a ■ v)} = 2u ■ v , (B5b) 

r,s 

holding for arbitrary vectors, u and v, with real-valued components. Together with the relations 

^ SrsXrCTi {3 ■ U) Xs = ^ ^rsXr (<? • «) (JiXs = Tr {(Jj (a • u)} = 2m (B6) 

r,s r,s 

and 

^ [Xr Xs] [Xr0"i Xs]* =Tr{(<7-u) {a-v) 2 a^ =2v 2 m, (B7a) 

XI [*> • ^) (<? • *0 Xs] [Xr (<? • u) OiXs]* = Tr | (a • u) 2 (<r • u) 2 o-jj = 2u 2 Vi , (B7b) 
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the photon self-energy for mixed time- and spacelike indices is evaluated to 



iU^(k,t,u) = 2e 5 



d 3 p j Pi + k 
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(B8b) 

Here it is important to note that in contrast to (|B6j) and ()B7j) . the index index i in (|B8aP and 
(|B8b|) is to be taken as covariant and not as Euclidean, i.e. {uj} = —u. Finally, we take into 
account that 

XI \Xr°i {° ■ U) Xs] [Xr<?j (ff ■ v) Xs]* = ^2 & r (° ' ^ aiXs ^ & r @ ' ^ 



r.s 



25ijU ■ v 



^2 [Xr(Ti (ff ■ U) Xs\ [Xr {ff ■ V) UjXs]* = ^ [Xr (<? ■ U) UiXs] [Xr<?j (<? " V) Xs]* 
r,s r,s 

= 2 (UiVj + UjVi — 5ijU ■ v) , 

which allows us to rewrite the photon self-energy for two space-like indices as 
iII<.(M,u) = 2e 2 |^3 {/y [(4. t (t)^(t)a^(«)/V, t (tt)+ 

^t(*) a p+^(*)^+fc4( u ) a P',t(«)_ 
- % [ a p-,t( t )/ 3 p+fc,^ t )^+^( n ) a P-t( n )+ 

^it(*)«i^(*)«jHE>)^it(«)4f] } • 

In order to keep the notation short, we have introduced 

£ = Pi(Pj + kj)+ pjjpj + h) + rjjjp- (p + k) 



(B9a) 
(B9b) 



(BIO) 



(Bll) 



p\p + k\ 

The verification of (|Bip in terms of wavefunction parameters is now straightforward. It follows 
immediately from (|36|) that 

d t (a^(t) ap+%A {t)) = ip^(t)a, +u {t) - i\p + k\a*^{t)^ +u {t) , (B12a) 
dt = jpo£ t (t)^(t) - i|p + k\^(t)a p+u (t) . (B12b) 
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Taking the time derivatives of (|B4p and (|B8ap and making use of relations (|B12ap and (|B12bp , we 
obtain 

d t iU^ (k,t,u) = ijfeMlI^A-.t.u) , (B13a) 

d t m<(k,t,u) = ikHUf^k^u) , (B13b) 

and hence 

d t pI&{k,t,u)-Ue i in&(iz,t,'u) = 0, (B14a) 

d t ill<(k, t, u) - ikHUf^k, t, u) = . (B14b) 

These are just the Ward-Takahashi identities (jBip . 



Appendix C: Remark on time dependence 



It has been shown in appendix [X] that keeping the correct sequence of the limits, first t — > oo 
and then e — > 0, is crucial for eliminating possible unphysical contributions from the vacuum 
polarization and that their interchange comes along with a divergent contribution from the latter. 
In this appendix, we consider the time evolution of ()137bp and (|137cp in order to highlight that 
keeping the correct order of limits is also essential to obtain physically sensible results from the 
mass-shift contribution (j!35bj) and to eliminate possible contributions from the interference term 
(I135cp . The latter is essential to obtain photon numbers that can be written as an absolute square 
and, as a consequence, are positive (semi-) definite. For this purpose, we first rewrite (|137bp and 
(fi~3TcD as 
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where we have introduced 

I^,(p,k,r,s,t) = 

J £ ^(p,k,r,s,t) = 

The time integral in (|C2bp is evaluated to 
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To handle (|C2aP for large times t > T S> r, we again take into account that the fermionic 
wavefunctions essentially have turned into (|129p in that range. Accordingly, we obtain 
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for t >T. To keep the notation short, we have introduced 
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If we now let e —> at some finite time, t > T, we see that ()C4j) turns into 
I £ ^(p, k, r, s, t) -> J M (£ fc, r, a, t) = / du j^ ST (p, fc, r, s, u)e iaJ £ M 
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with g(u!,T,t) given by 
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>(t-T) 



On the other hand, (|C3p does not vanish, but behaves like 

J»(P, k , r , s, t) = eu c {p, r)-f^v c {p + k, s) — + 0(e) , for e -)• . 
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Hence, when first taking the limit e — > at some finite time, t, expression (jClbp does not vanish. 
Furthermore, we can infer from Fig. [27] that it generally persists in the limit t — > oo. 

If one first takes t — > oo, however, it follows from (|C3P and (|C4p that these expressions turn 
into (|14U|) and (|143|) . respectively, so that (jClbp vanishes when taking the subsequent limit e — > 0. 

In addition to the elimination of possible unphysical contributions from (|135ap and (|135cp . 
keeping the correct sequence of limits is also crucial to obtain physically reasonable results for the 
contribution from ()135bp describing mass-shift effects only. To illustrate this, we show the time 
evolution of (jClap for different values of e in Fig. 
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FIG. 27: Time evolution of the interference contribution (jClbft for oj£ = 1 GeV and different values of Ac*. 
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FIG. 28: Time evolution of the mass-shift contribution (IF2[) for different values of e for m2{t) (left panel) 
and 7713(f) (right panel). 



Taking first e — > at some finite time, i, corresponds to the curve labeled by e = in both 
panels. If we then take t — > 00, we obtain an asymptotic value for (jClap . which is by several 
orders of magnitude larger than the one against which the asymptotic values of (|Cla|) for finite e 
converge in the limit e — > and which accordingly corresponds to (|62|) . Hence, we have seen that 
interchanging both limits also leads to unphysical results from the pure mass-shift contribution. 
Moreover, we see that for sufficiently small values of e, expression (|Clap for finite times t > T ^> r 
coincides with its asymptotic value for e = 0. Consequently, the interpretation of (jClaj) as a photon 
number at finite times is doubtful as well. 

In order to highlight these aspects in particular, we take a closer look at (|C6p and (|144p . 
We see that the integral expressions coincide in both cases whereas the frequency expressions 
entering (|C6p feature an additional term of e luJit jiuji (i = 0...4), which has been eliminated in 
the corresponding frequency expressions entering ()144p after switching off the electromagnetic 
interaction adiabatically. 

In order to see why this additional term is artificial in each case we again consider ()ClaP for the 
correct sequence of limits, but instead perform an instantaneous switching of the electromagnetic 
interaction at t = =Fl/e, i.e., 

Hem (i) -> 6 e (t)H EM (t) , with 6 £ {t) = 6 ( \ - t 2 ) and e > . (C9) 
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In this case, (|143p is replaced by 
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In order to keep the notation more shorthand, we have defined 



T 



-AujT 



iuj 



(CIO) 



(Cll) 



It follows from ([C7]) and (jCTTj) that taking the limit e -> in (fCTO]) has the same effect as taking 
t — > 00 in (|C6p . This can also be inferred from Fig. 1291 showing the time evolution of (jClap for 
different values of e with the electromagnetic interaction being switched according to (|C9p , 
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FIG. 29: Time evolution of the mass-shift contribution ([Clap for 777,2(7;) (left panel) and 7713(7;) (right panel) 
with the electromagnetic interaction being switched in terms of (|C9[) . If we take the successive limits t — > 00 
and then e — > 0, we obtain the same result for (|Cla[) as if we had kept the regularization P%|) but interchanged 
both limits instead. The corresponding value is indicated by the dashed horizontal line in each case. 



Thus, if we switch the electromagnetic interaction according to (|44p but take the incorrect 
sequence of limits we obtain the same result for (jClaP as if we took the correct sequence of limits 
but switched the electromagnetic interaction according to (|C9p . This shows us that the former 
case is formally equivalent to keeping the exact sequence of limits for an unphysical instantaneous 
switching (|C9p . 

In order to emphasize why a switching of the electromagnetic interaction by means of (|C9p is 
unphysical, we consider again the Ward-Takahashi identities (|6ip for the photon self-energy. We 
have shown in section [TT] and also demonstrated in greater detail in appendix[B]that these identities 
are fulfilled when coupling the quark field to a time-dependent scalar background field, (j>(t). But 
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one has to keep in mind that a switching of the electromagnetic interaction leads to a corresponding 
modification of the photon self-energy, i.e., 

in< (k,h,t 2 ) in^f (M!,t 2 ) = f e (t 1 )f e (t 2 )iu< u (k,t 1 ,t 2 ) , (C12) 

with the electromagnetic interaction then being fully persistent at any given time, t. Because of 
the additional factor of f e {ti), relation ([61 p is not fulfilled in the first place, but we instead have 

d tl m<> E (k,h,t 2 ) - i^'in<f (M l5 i 2 ) = -f E (ti)f e (t2)w££$,ti,t 2 ) . (Ci3) 

However, if the electromagnetic interaction is switched according to (|44p . the time derivative reads 

f e (t) = e sign(t)e- £ l*l , (C14) 

and, as a consequence, the r.h.s. of (|C13p vanishes for e — > such that ([6T]) is restored again. 
To the contrary, if we switch the electromagnetic interaction by means of (|C9p . we have no such 
restoration since 

£ (t) = 6 ( - - t J 5 ( - + 1 J - ( - +t J 5 ( - - t J . (C15) 



Hence, switching the electromagnetic interaction according to (|C9p effectively leads to a violation 
of the Ward-Takahashi identities. Since (|C9P leads to the same asymptotic photon numbers as 
(|44|) with the limits being interchanged, we see that the latter procedure also leads to an effective 
violation of (EH). 



Appendix D: Determination of Bogolyubov parameters £p, s (t) and r]p s (t) 



In this appendix, we show explicitly that the system of equations (|84|) for the Bogolyubov 
parameters, £p S (t) and rjff s (t), is solved by (|85p . For this purpose, we first rewrite the fermionic 
field operator, ip(x), in terms of the transformed operators (|71|) . i.e., 



d 3 p 
(2^)3 



(Dl) 



Here we have made use of ([76]) and introduced the transformed positive and negative energy 
wavefunctions 



i>p,s,i(x) = £fr s (t)^p,si{x) - r)p,s{t)i)$rf(x) . 
It follows from ()72p that the inverse transformation relations are given by 



(D2a) 
(D2b) 

(D3a) 
(D3b) 



Now the strategy is to rewrite ()84|) as matrix-element conditions for ipp tS ti(x) from which the latter 
are determined. The corresponding Bogolyubov parameters, £,p s (t) an< ^ Vpsit), are then given by 

CpA f ) = ^1,sA x )$p,sA x ) = v^WVy,^) = ^p, s ^ x )^pM x ) > ( D4a ) 

Vfrsit) = -^p,st( x )$p,sl( x ) = ^sA^pM 00 ) = $p,st( x )^P,sl( x ) ■ ( D4b ) 
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Expressing (|68ap - (|68b|) in terms of (|D3ap - (jD3bp and making use of (|72p . we can rewrite (|84p in 
terms of ipp,sfi{^) as 



To keep the notation short, we have introduced 

hD(t) = -i7o7 • V + 7om(t) . 

Taking into account that 

which follows from (|D2[> and (|68ap . one can show that ([D5aJ)-(jD5bJ) are equivalent to 

V'^ s ,t( a; )^- D ( t )V'p>s,t( a; ) = . 







(D5a) 
(D5b) 

(D6) 

(D7) 

(D8a) 
(D8b) 



The Bogolyubov transformation (|7ip hence corresponds to a reexpansion of the fermion-field oper- 
ators in terms of the instantaneous eigenstates of the Hamilton-density operator (|D6p . In order to 
solve Eqs. (|D8ap - ()D8bp for the matrix elements, we parametrize Vy,st4-( x ) analogously to (f35j) as 



7 ( \ - ( %.t4.(*)x« 



(D9) 



Writing (|D2aP and (|D2bp in terms of wavefunction parameters, it follows from (|38p and (|72p that 



Furthermore, relations 



imply 



1 . 



(D10) 

(Dlla) 
(Dllb) 



Upon insertion of (|D9[) into (|D8ap - ()D8bp and making use of (|Dllap -( )Dllbp the diagonalization 
conditions can be written as 



P (c^, t (*)/%;t(t) + %it(*)%(*)) + m ( t ) (l«P,t(*)| " 

p (^| t (t) - a| t (t)) + 2m(t)d p - it (t)4- jt (t) 
Expressing the wavefunction parameters in terms of polar coordinates, i.e., 

a p - )t (i) = la^^le 1 ^ , 



%(t) , 



(D12a) 
(D12b) 

(D13a) 
(D13b) 



and taking into account that the l.h.s.s of both (|D12ap and (|D12bp are real valued, one obtains 

smiMt) -<f>p(t)] = , 
cos[0 Q (t)-<^(i)] = (-1)" . 
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Therefore, we can rewrite (|D12aP - (|D12b|) as 



(-l) n 2p\a^(t)\ % t (i) +m(t)[\a p>t (t)\ - (5^{t) 



(-l) n pM5ftt(*)| - %,t(*) J +2m(t) |«p, t (t)| /%, t (t) 
which is equivalent to 



, 



l«at(*)l ~ /%,t(*) 



m(t) 
E^j 



(D14a) 
(D14b) 

(D15a) 
(D15b) 



Relation (|D15aP implies that n must be even. Thus, without loss of generality, we can assume 
that cp a (t) = (f>s(t) = <j){t). Together with the normalization condition (|D10j) . it follows from 
(|D15ap - (|D15bp that we have 



«#,t(*) 



'Ep(t) + m(t) ; 



2EJt) 



< Ep(t)-m(t) j m 
2E p (t) 



(D16a) 
(D16b) 



Since ipp s ^A x ) an d ^p,s,t4.(' E ) coincide with each other for t < t , <f)(t) can, in principle, be chosen 
freely as long as it fulfills the asymptotic condition (j){t) — > E^t for t —> — oo. It is hence convenient 
to choose <f)(t) = Ep(t)t so that our Bogolyubov transformation corresponds to an expansion of 
ip{x) in terms of positive and negative energy wavefunctions of the respective current mass, m(t). 
As a consequence, we finally have for the transformed positive and negative energy wavefunction 
parameters 



5?u(*) 



< ' Ep(t) + m(t) gi%(t)f 
2Ejf(t) 



l Epjt) -m(t) iEf(t)t 
2Ep(t) 



<Ep(t) - m(t) „_ i%(t)t 



2E p {t) 



(t)t 



2E p {t) 



(D17a) 
(D17b) 
(D17c) 
(D17d) 



and hence 



,\E p {t)t 



1 Ep(t) + m{t) 
2Ep(t) 

*EJt) + m(t) 



2E p (t) 

These are just the Bogolyubov parameters 



<y,f(*) + 



a pA f ) + 



lEp{t)-m(t) 
2Ep(t) 

1 'Epjt) - m(t) 
2EJt) 



PpAt) 



PpAt) 



(D18a) 
(D18b) 
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Appendix E: Evaluation of k) 

In this appendix, we demonstrate how to evaluate Iij(p,k) = 1..4), which characterize the 
direct contributions from first-order QED processes (equal indices) and their interference among 
each other (different indices). We will carry out the calculations explicitly for one direct contribu- 
tion and one interference contribution, respectively. The evaluation of the remaining contributions 
proceeds analogously. In order to calculate the direct contribution from the spontaneous creation 
of a quark-antiquark pair together with a photon, In(p,k), we first recall that the corresponding 
transition amplitude is given by 

k, r, s) = ieapa- %u b (p, r)^v b (p + k,s)( — — -=- J . (El) 

\w?(p,fc) (p, k) J 

In order to carry out the spin summation in In(p, k), we take into account that 

2^ u b (p, s)u b (p, s) = 2 Eb . ( E2a ) 

s V 

Y,v b {p,s)v b {p,s) = , (E2b) 



2E} 



where p^ = (E^, —p). Hence, In(p, k) can be rewritten as 



p+k V 



where we have introduced = {E\ ?,p + k). Using standard Dirac trace techniques, we find 
that 

{lu($ ~ m b )^(f + m b )} =4:[q l/ p ll + q fl p u -r] l/fl (q-p + ml)} . (E4) 

Taking into account that ^"(k) has a purely spacelike structure given by (|48p . carrying out the 
contraction entering (|E3p finally yields 



,-. r , 99 o pxipx + cjjA + mi \ f 1 1 \ . . 

k f) " 2e ^ [ + % + k ) " w$ ) ' <E5) 

The evaluation of the direct contributions from quark and antiquark bremsstrahlung and quark- 
antiquark annihilation into a photon is performed in the same way. Next, we turn to the evaluation 
of the interference contributions. In its course, we encounter spin summations over mixed tensor 
products between positive- and negative-energy spinors. These evaluate to 

2^ u b {p, s)v b (p, s) = — b , (E6a) 



v 



v b (p, s)u b (p, s) = - ^ p<y \ * ^ 



24 

where we have introduced rut = (0,p/p) and taken into account that 



(E6b) 



v 

ipu b (p,s) = -v b (p,s) , (E7a) 
ftpv b (p,s) = u b (p,s) . (E7b) 
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With help of these relations, the contribution describing the interference between the spontaneous 
creation of a quark/antiquark pair together with a photon and quark pair annihilation into a 
photon, Ii/i(p, k,t), can be rewritten as 

k) =e . . — ^ - — — =- — _ + 



4E ? Eb p+k \ W l(P>^) W l(P:*0/ W(£fc) ^l(p,k)J 

• Tr { 7l/ (j* + £ + mMfrtffM + m 6 )jt f } (E8) 

Thus, in contrast to the direct contributions, carrying out the spin summations for the interfer- 
ence contributions gives rise to traces over products of six Dirac matrices. Making use of the 
anticommutation relation 

and cyclic trace invariance, these traces can be reduced to 

Tr{ 7 Q 7 /3 7V7 e 7 f } = ^Tr ^VtV} ~ < 7 Tr {tVtV} + 

rf 5 Tr {7V7V} - ^'Tr {7V7V} + 
7? °CTt{ 7 ^7 7 7V} • (E9) 
With the help of (|E9[) . the trace entering (|E8[> is evaluated to 
Tr {7^ + ^ + m^i^^f + m b )^} 

= - — ^ pi(pj + fcj) + pj(pi + fej) + rjijp- (p + k) + r/jjp|p + k\ . 
p\p + k\ 1 i 



(E10) 



We have considered only spacelike indices since the trace expression entering (|E8j) is still contracted 
with 7 / " / (fc) which has an exclusively spacelike structure. This contraction leads to 



7 ^(fc)Tr {7^ + jg + m 6 )# ?+ ^(^ + ™{-¥ p -} 

|p + fc| 



V + *T 1 + — ^vf^ I . (EH) 




so that (jE8J) finally reads 

hM k)=~ 2e 2 " "J: P+k P \p + k\ 1 + 



|i?+£| s 



p p+fc 
1 1 




(E12) 



+ 



u)\(p,k) u)%(p,k)J \ul(p,k) uf(p,k) 



The evaluation of the remaining interference contributions follows the same steps. After evaluating 
all direct and interference contributions the substitutions p — )• p — k and x — > — 1 show that the 
remaining integration over d 3 p yields the same contribution for l22(p,k,t) and Iss(p,k,t), for 
Ii2(p,k,t) and Iiz(p,k,t) and for l2i(p,k,t) and l3<i(p,k,t). For that reason, theses contributions 
have been taken together as shown in (|128|) in each case. 
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Appendix F: Remarks on mj = for an instantaneous mass shift 



In this appendix, we discuss the special case of mb = for an instantaneous mass shift. It re- 
quires some special considerations since the frequencies describing quark/antiquark bremsstrahlung 
and pair annihilation into a photon are no longer negative or positive definite. We first take into 
account that we have 



ap = -j= (cos (pjf + sin p£) 
Pp = -j= ( cos <Pp- sin *Pp) 



(Fla) 
(Fib) 



for rrib = 0. The asymptotic direct contributions hence turn into: 



hi(p,k) = — 



h2(p,k) 



1 + 



\p + k\ 

E c - 
p+k . 

1 



1+ 4 



i + 
i 



x(px + Ur 
\p + k\ 

\ 2 



\p + k\ + p + uj: 
, \P + k\ 



p+k P fc, 



p+k 
1 



\p+k\ -p-uj^ 
, \P + k\ S 



1 + 



+ 



p_ 



i 



x{pX + OJl 

+ 
\ 2 



p+fc P k. 



E c 



1 



p+fc 
1 



P/ 



+ 



1 + 
1 



x(px +LO^) 

W+k\ 

\ 2 



(F2a) 



(F2b) 



(F2c) 



Furthermore, it follows from (|127|) that I\2{Pi k) and l24,(p, k) vanish for m& = and the remaining 
interference contributions simplify to 



Iu(p,k) 



h3(p,k) 



2E c :E c - 

P p+k 



2E c :e° 

P p+k 



1 + 



x(px + ojj:) 
W+k\ 



+ 



1 



Jp+ k\ + p + u f : 
1 



p+k P k, 



1 



\p+ k\ - p-u- k 
1 



+ 



E c -, + E c , + u t 

p+k P fc, 



\p + k\ -p + ujj: . . 



(F3a) 



(F3b) 



Since the frequencies describing quark/antiquark bremsstrahlung and pair annihilation into a pho- 
ton are no longer negative or positive definite, we have to investigate, whether the loop integral 
features possible infrared {p — > 0) and/or (anti-) collinear (x —> ±1) singularities. In the limit 
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p — > 0, they behave as 

e 2 (l + x) 



hi(p,k) 
Iu(p,k) 




1 + 



UJt 



uji + m 2 




ul + m 2 c + EZ + Ut 



1 



ujl + m 2 ) P 2 (l-*) 



1 



col + m 2 / P 2 (l + *) 



oil 



2 P A /w| + m2 




+ 0(p) , (F4a) 
(F4b) 
(F4c) 
(F4d) 
(F4e) 



Since the integration measure d 3 p contributes an additional factor of p 2 to each contribution, no 
infrared singularities will show up in the loop integral for mt, = 0. As a next step, we consider the 
collinear limit x — > 1. In this limit, the direct contribution from quark/antiquark bremsstrahlung 
turns into 



W+k\ 
E c - 

p+k . 



x=l 



1 + 

P 



p \ (p + ujj:) 2 +P 2 



p 2 uj"l(x — 1) 



+ 0(1) 



(F5) 



with all other contributions staying finite. Furthermore, the direct contribution from pair annihi- 
lation into a photon features an anti-collinear singularity in the limit x — > 1 if p < wj, i.e., 



Iu(p,k) -> — 1 



|p + k\ 



E< 



p+k 



(F6) 



For p > wr, i44(p, fe, t) stays finite for x — )• — 1. All other contributions generally stay finite in that 
limit. Therefore, in the limit m& — )■ 0, the loop integral develops a collinear and an anticollinear 
singularity for the contributions from quark/antiquark Bremsstrahlung and quark pair annihilation 
into a photon, respectively. 
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